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We review briefly the quantum fidelity approach to quantum phase transitions in a pedagogical 
manner. We try to relate all established but scattered results on the leading term of the fidelity into 
a systematic theoretical framework, which might provide an alternative paradigm for understanding 
quantum critical phenomena. The definition of the fidelity and the scaling behavior of its leading 
term, as well as their explicit applications to the one-dimensional transverse-field Ising model and 
the Lipkin-Meshkov-Glick model, are introduced at the graduate-student level. In addition, we 
survey also other types of fidelity approach, such as the fidelity per site, reduced fidelity, thermal- 
state fidelity, operator fidelity, etc; as well as relevant works on the fidelity approach to quantum 
phase transitions occurring in various many-body systems. 

PACS numbers: 03.67.-a, 64.60.-i, 05.30.Pr, 75.10.Jm 



Contents 



E. Density-functional fidelity 



I. Introduction 

A. Overview: quantum phase transitions 

B. Brief historical retrospect 

C. About the review 



il 
i 



II. Quantum fidelity: a measure of similarity 
between states 4 

A. Pure state and mixed state fidelity |4 

B. Quantum state overlap and adiabatic 
evolution 

III. Fidelity and quantum phase transitions 

A. Quantum phase transitions: fidelity 
perspective 

B. Example: the one-dimensional 
transverse-field Ising model 



IV. 



Fidelity susceptibility, scaling and 
universality class [ll| 

A. The leading term of the fidelity and dynamic 



structure factor 

B. Scaling analysis and universality class 

C. Example A: the one-dimensional 
transverse-field Ising model 

D. Example B: the Lipkin-Meshkov-Glick 
model 

E. Higher order of the fidelity 



V. 



Fidelity per site, mixed state fidelity, and 
related 2£ 

A. Fidelity per site [Ic 

B. Reduced fidelity 

C. Fidelity between thermal states 

D. Operator fidelity 



'Electronic address: ' sjgu@phy.cuhk.edu. hk| 



VI. Fidelity in strongly correlation systems |24 

A. Pure-state fidelity 24 

1. One-dimensional spin systems 24 

2. One-dimensional fermionic systems |27 

3. The fidelity in topological quantum phase 
transitions 

B. Mixed-state fidelity 

1. Reduced fidelity in quantum phase 
transitions [35 

2. Thermal state fidelity in strongly 
correlated systems [SS 

VII. Numerical methods for the ground-state 

fidelity [si 

A. Exact diagonalization 34 

B. Density matrix renormalization group |3f 

VIII. Summary and outlook [36 

Acknowledgments [36 

References [H 



I. INTRODUCTION 

A. Overview: quantum phase transitions 

Quantum phase transitions [l| of a quantum many- 
body system are characterized by the change in the 
ground-state properties caused by modifications in the 
interactions among the system's constituents. Contrary 
to thermal phase transitions where the temperature plays 
a crucial role, quantum phase transitions are completely 
driven by quantum fluctuations and are incarnated via 
the non-analytic behavior of the ground-state properties 
as the system's Hamiltonian H{X) varies across a transi- 
tion point Ac. 
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From the point view of eigenenergy, quantum phase 
transitions are caused by the reconstruction of the Hamil- 
tonian's energy spectra, especially of the low-lying exci- 
tation spectra More precisely, the low-energy spec- 
tra can be reconstructed in two qualitatively different 
ways around the critical point Ac , and hence the physical 
quantities show different behaviors. The first one is the 
ground-state level-crossing in which the first derivative 
of the ground-state energy with respect to A is usually 
discontinuous at the transition point. Such a transition 
is called the first-order phase transition. The second one 
corresponds roughly to all other cases in the absence of 
the ground-state level-crossing. It is usually a continuous 
phase transition. 

Traditionally, continuous phase transitions can be 
characterized by the Landau-Ginzburg- Wilson sponta- 
neous symmetry-breaking theory where the correlation 
function of local order parameters plays a crucial role. 
Nevertheless, some systems cannot be described in this 
framework built on the local order parameter. This 
might be due to the absence of preexistent symmetry in 
the Hamiltonian, such as systems undergoing topological 
phase transitions [a| and Beresinskii-Kosterlitz-Thouless 
phase transitions [¥, ^ . 



B. Brief historical retrospect 

In recent years, ambitions on quantum computer 
and other quantum information devices have driven 
many people to develop quantum information theory Q . 
Though a practicable quantum computer seems still a 
dream, progresses in quantum information theory have 
developed other related fields forward. A noticeable one 
is the relation between quantum entanglement and quan- 
tum phase transitions 0, i, i, [13, El, [13, [13, M, [H, 
[igI . [itI [TsI . [T9| . Since the entanglement is regarded as 
a purely quantum correlation and is absent in classical 
systems, people think that the entanglement should play 
an important role in quantum phase transitions. Though 
a unified theory on the role of entanglement in quantum 
phase transitions is still unavailable, some definitive con- 
clusions have been commonly accepted 20]. 

Another attractive approach is the quantum fid elity 

[U, M M MM MM M M M M, M M 

a concept also emerging in quantum information the- 
ory. The fidelity measures the similarity between two 
states, while quantum phase transitions are intuitively 
accompanied by an abrupt change in the structure of the 
ground-state wavefunction, this primary observation mo- 
tivates people to explore the role of fidelity in quantum 
phase transitions jssl . [s^ . Since the fidelity is purely a 
quantum information concept, where no a priori knowl- 
edge of any order parameter and changes of symmetry of 
the system is assumed, it would be a great advantage if 
one can use it to characterize the quantum phase tran- 
sitions. Many works have been done along this stream 

M M M m M M El, Si, S, H S, M M, M Si, 
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The motivation of the fidelity approach to quantum 
phase transitions can be traced back to the work of 
Quan et al [35] in determining two ground-state phases 
of the one-dimensional transverse-field Ising model by 
the Loschmidt echo. The Loschmidt echo (IIS] has been 
introduced to describe the hypersensitivity of the time 
evolution to perturbations experienced by the environ- 
mental system. They found the quantum critical be- 
havior of the environmental system strongly affects its 
capability of enhancing the decay of Loschmidt echo. 
Since the Loschmidt echo is defined as the overlap be- 
tween two time-dependent states corresponding to two 
points separated slightly by a target spin with Ising in- 
teraction, its decay around the critical point represents 
a large distance between two states. Subsequently, Za- 
nardi and Paunkovic (36| proposed out that a static fi- 
delity might be a good indicator for quantum phase tran- 
sitions with examples of the one-dimensional transverse- 
field XY model and the Dicke model. Similar idea was 
also proposed by Zhou and Barjaktarevic [13]. Motivated 
by these works, the fidelity approach to quantum phase 
transitions was quickly applied to free fermionic systems 
[3^ and graphs f39'] , matrix-product states [4^ , and the 
Bose-Hubbard model ^]. An attempt to understand 
quantum phase transitions from the thermal fidelity was 
also made (4^. At that time, the successes of the fi- 
delity in these studies [H, M M [M, l39, 4l, M [H 
gave peoples a deep impression that the fidelity is able 
to characterize any quantum phase transition, includ- 
ing those cannot be described in the framework of 
Landau-Ginzburg- Wilson theory, such as the Beresinskii- 
Kosterlitz-Thouless transtions and topological transi- 
tions. Two groups addressed the role of the leading term 
of the fidelity in the quantum critical phenomena. Za- 
nardi et al introduced, based on the differential-geometry 
approach, the Riemannian metric tensor |^44] inherited 
from the parameter space to denote the leading term in 
the fidelity, and argued that the singularity of this metric 
corresponds to quantum phase transitions. While You et 
al introduced another concept, the so-called fidelity sus- 
ceptibility (FS) [43], and established a general relation 
between the leading term of the fidelity and the struc- 
ture factor (correlation functions) of the driving term in 
the Hamiltonian. Both of them obtained also that, if one 
extend the fidelity to thermal states, the leading term 
of the fidelity between two neighboring thermal states is 
simply the specific heat. In the following, we will use 
"fidelity susceptibility" to name the leading term of the 
fidelity because it not only denotes mathematically the 
fluctuation of the driving term, such as the specific heat 
derived from the internal energy, but also is closer to the 
picture of condensed matter physics, i.e. the response of 
the fidelity to driving parameter. From then on, the field 
of the fidelity approach to quantum (or thermal) phase 
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transitions can be divided roughly into two streams. The 
first stream still focuses on the fidelity itself, for which 
the distance between two points in the parameter space is 
still important, while the second stream pays particular 
attention to the leading term of the fidelity. 

Along the first stream, a connection between the fi- 
delity, scaling and renormalization was introduced by 
Zhou ^45. . .46.] , in which the fidelity between two reduced 
states of a part of the system described by a reduced- 
density matrix was proposed. Zhou et al [s^ tried to 
understand the fidelity from a geometric perspective. In 
works of Zhou and his colleagues, the fidelity is aver- 
aged over the system size, and is named as fidelity per 
site. They found that the fidelity per site is a very use- 
ful tool for various interacting systems. Interestingly, 
the fidelity per site, as an analog of the free energy per 
site, can be c omputed in the context of tensor network 
algorithms [Ei, HS, [l43] . 

While along the second stream, several questions ap- 
peared at that time. 1) Since the leading term of 
the fidelity is a combination of correlation functions, 
which seems a tool widely used only in the Landau- 
Ginzburg- Wilson theory, is the fidelity still able to de- 
scribe the Beresinskii-Kosterlitz-Thouless and topologi- 
cal phase transitions? 2) What is the scaling behavior of 
the fidelity and its relation to the universality class? 3) 
How about the thermal phase transitions and those quan- 
tum phase transitions induced by the continuous ground- 
state level-crossing where the perturbation method is not 
applicable. Most subsequent works are more or less re- 
lated to these questions, though some topics are still con- 
troversial. 

Based on the general relation between the leading term 
of the fidelity and correlation functions of the driving 
term , Venuti and Zanardi [J?! applied the traditional 
scaling transformation, and obtained an interesting scal- 
ing relation between the dynamic exponent, the dimen- 
sion of the system, and the size exponents of the fidelity. 
A similar scaling relation was also obtained numerically 
by Gu et al 49] in their studies on the one-dimensional 
asymmetric Hubbard model. Both relations imply that 
the fidelity susceptibility might not have singular be- 
havior in some cases, such the Beresinskii-Kosterlitz- 
Thouless transition occurring in the asymmetric Hub- 
bard model at half-filling [4§] . 



On the other hand, Yang [54] tried to understand 
the singular behavior of the fidelity susceptibility from 
the ground-state energy density and pointed out that 
the fidelity susceptibility might not be able to detect 
the high-order phase transitions. A little surprising is 
that their example, i.e. the effective model of the one- 
dimensional XXZ chain, which undergoes a Beresinskii- 
Kosterlitz-Thouless transition of infinite order at the 
isotropic point, shows singular behavior in the fidelity 
susceptibility. Similar analysis on the Luttinger Liquid 
model with a wave functional approach was also done by 
Fjrestad [6(|. The further investigations on spin-1 XXZ 
chain with uniaxial anisotropy by Yang et al [ssl . pTOj sup- 



ported partially their previous conclusion and the scaling 
relation obtained by Venuti and Zanardi ^47i] . 

Later, Chen et al addressed the feasibility of the fi- 
delity susceptibility in quantum phase transitions of vari- 
ous order by the perturbation theory, and concluded that 
the fidelity susceptibility cannot describe the phase tran- 
sition of infinite order. This conclusion conflicts with 
both Yang's works on the one-dimensional XXZ model 
[5^ and the subsequent studies on the one-dimensional 
Hubbard model [63|, but supports previous conclusion 
obtained by You et al^^. Therefore, the issue on fi- 
delity in describing high-order phase transitions seems 
still controversial. 

Recently it was realized that the fidelity susceptibility 
can be either intensive, extensive, or superextensive, then 
the critical exponents of the rescaled fidelity susceptibil- 
ity at both sides of the critical point can be different [t^ . 
In addition to the fidelity susceptibility, the sub-leading 
term of the fidelity might appear when parameters are 
changed along a critical manifold ff^ . 

It became a branch of the story when Hamma et al 
65] firstly touched the feasibility of the fidelity in topo- 
logical phase transitions. They found that though the 
fidelity shows an obvious drop around the critical point 
of a topological transition, it cannot tell the type of tran- 
sition. Almost one year later, three groups revisited the 
role of fidelity in the topological transitions. Zhou et al 
[m studied the fidelity in the Kitacv honeycomb model 
and found that fidelity has shows singular behavior at the 
critical point. Yang et al [63] studied the fidelity suscep- 
tibility in the same model and obtained various critical 
exponents, they also witnessed a kind of long-range cor- 
relation in the ground state of Kitaev honeycomb model. 
While Abasto et al ^68:] studied the fidelity in the de- 
formed Kitaev toric model and obtained a form of fidelity 
between thermal states. The three groups drew a similar 
conclusion that the fidelity can describe the topological 
phase transitions occurring in the both models. 

A noteworthy advance in the fidelity approach is the 
success of using the state overlap to detect quantum 
critical point by a nuclear-magnetic-resonance quantum 
simulator (tJ. It was observed that the different types 
of quantum phase transitions in the transverse-field Ising 
model can be witnessed in experiments. Such an advance 
is remarkable. It makes the fidelity approach to quantum 
phase transitions no longer purely theoretical. 

On the other hand, the global-state fidelity cannot 
characterize those quantum phase transitions induced by 
continuous level-crossing due to its collapse at each cross- 
ing point. Kwok et al ^] firstly tackled this type of 
phase transition with the strategy of the reduced fidelity, 
which actually was introduced in previous works [45l [Sa ] . 
Meanwhile, Ma et al [TSJ also studied the critical behav- 
ior of the reduced fidelity in the Lipkin-Meshkov-Glick 
model. The reduced fidelity was latter applied to the 
one-dimensional transverse-field Ising [^, and XY 
models !)() , the dimerized Heisenberg chain 82\ , and the 
one-dimensional extended Hubbard model |88i] . 
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Despite of the absence of the thermal phase transition 
in the one-dimensional XY model, the thermal-state fi- 
delity was firstly used to study the crossover occurring 
in the low-temperature critical region [4^ . Interestingly, 
the leading term of the thermal-state fidelity was later 
found to be just the specific heat [H, The thermal- 
state fidelity was also applied to the BCS superconduc- 
tivity and the Stoner-Hubbard model with the mean-field 
approach [s^. Moreover, Quan and Cucchietti [zll tried 
to find the advantages and disvantages of the fidelity ap- 
proach to the thermal phase transitions. 

Finally, though we focus on the fidelity between the 
static ground state only, we would like to mention that 
the Loschmidt echo has also been widely applied to study 
the quantum phas e transitions [sl, H M, M, 113, H M, 
[Tonl . flQli flQl - Iloi [foi . [TnHj . in studies of the Loschmidt 
echo, one needs to consider the dynamic behavior of the 
fidelity of a target object, for instance, a spin coupled 
with all other spins in the Ising chain. Then the deco- 
herence property should be taken into account. These 
issues are beyond the scope of this review. 



C. About the review 

The main purpose of this review is to gather these 
distributed works into a unified paradigm, then pro- 
vides interested readers, especially beginners, a system- 
atic framework of the fidelity approach to quantum phase 
transitions. Some practical and numerical methods, such 
as the exact diagonalization and density matrix renor- 
malization group, are introduced too. We try to keep 
the treatment as simple as the subject allows, showing 
most calculations in explicit detail. Since the field is still 
quickly developing, such a review is far from complete- 
ness. We hope that the article can offer some introduc- 
tory essays first, then to arouse more wonderful ideas. 

The article is organized as follows. In Section |TT1 we 
give a brief overview on the fidelity measure and its prop- 
erties in an adiabatic evolution exampled by a 1/2 spin 
subjected to an external field. In Section IIIIl we intro- 
duce in considerable detail the general relations between 
the fidelity and quantum phase transitions, and try to 
illustrate the role of fidelity in quantum phase transi- 
tions by the one-dimensional transverse- field Ising model. 
In Section IIVI we focus on the leading term of the fi- 
delity, i.e. the fidelity susceptibility, and discuss its gen- 
eral properties around the critical point. We also use 
the one-dimensional transverse-field Ising model and the 
Lipkin-Meshkov-Glick model as examples for the fidelity 
susceptibility in describing the universality class. In Sec- 
tion El we review other types of fidelity in the quantum 
phase transition, such as the fidelity per site, partial- 
state fidelity, thermal-state fidelity, operator fidelity, and 
density-functional fidelity. In Section IVIi we give a sur- 
vey on the fidelity approach to quantum phase transitions 
in various strongly correlated system. In Section [VI II we 
show how to calculate the fidelity and fidelity suscepti- 



bility via some numerical methods. An outlook and a 
summary will be presented in the concluding section. 

II. QUANTUM FIDELITY: A MEASURE OF 
SIMILARITY BETWEEN STATES 

In this section, we introduce briefly the concept of 
quantum fidelity and discuss its properties in a simple 
quantum-state adiabatic evolution of a 1/2 spin sub- 
jected to an external field. 

A. Pure state and mixed state fidelity 

In quantum physics, an overlap between two quantum 
states usually denotes the transition amplitude from one 
state to the another [Hi HHHHll. While from the 
point view of information theory, the overlap can measure 
the similarity (closeness) between two states [1^ [2^ [13] . 
That is the overlap gives unity if two states are exactly 
the same, while zero if they are orthogonal. Such an in- 
terpretation has a special meaning in quantum informa- 
tion theory [al sinc e physicists in the field (for examples, 
Ref |l06l . IIOTI . [10_8|) hope that a quantum state can be 
transferred over a long distance without loss of any in- 
formation. The overlap between the input and output 
states becomes a useful measure of the loss of informa- 
tion during the transportation. The overlap is used to 
define the fidelity in quantum information theory. 

To be precise, if we define the overlap between two 
pure states as 

/(vl/',vl/) = (vlr'|vl,), (1) 

the fidelity is simply the modulus of the overlap, i.e. 

F(vl/',vl/) = |(vl/'|vl/)| (2) 

where j^*), j^*') are the input and output states respec- 
tively, and both of them are normalized. The fidelity has 
a geometric meaning as well. Since a pure state in quan- 
tum mechanics mathematically is a vector in the Hilbert 
space, then according to Linear algebra, an inner product 
of two vectors a, b is 

a • b = a6cos(6') (3) 

where a{b) is the magnitude of a(b), and 9 is the angle 
between them. In quantum mechanics, wave functions 
are usually normalized, and the fidelity represents the 
angle distance between two states. 

The fidelity has the following expected properties {ax- 
ioms) ^ 

< F(*',^') < I, (4) 
F(^'',*) = F(*,«''), (5) 
F(C/^'',C/*) = F(*',«'), (6) 
F(*i®*2,*'i®*2) = F(4''i,*i)F(*^,*2), (7) 
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where U denotes a unitary transformation and ^'1(2) is 
the state of one subsystem. For pure states, the global 
phase difference may affect the overlap, but not the fi- 
delity. 



Example: 

The quantum state of a single spin can be ex- 
pressed in the basis {| t), I i)}- For two nor- 
malized states of the spin, say 

1^(0)) = cos^l t) +sin0| i), 
\^{6')) = cos^'l T)+sin0'| I), 

the fidelity between them is 

F{^{e'),^{0)) = \ cos{e~0')\. 



The quantum fidelity between two mixed states (p, p') 
is defined as UM 



F{p, p') = trVpi^W^ 



(8) 



Here p{p') is semi-positive defined and normalized, i.e. 
trp =tTp' = 1. The definition satisfies the expected prop- 
erties of the fidehty, i.e. Eqs. (gllTl). 

It is not easy to evaluated the fidelity between two arbi- 
trary mixed states. Nevertheless, there are some special 
useful cases: 

1) If both states are pure F{p,p') = 

2) If one of state is pure, i.e. p = |\I')(\I'|, then 
F{p, p') — ^/{^\f7\W}, which is simply the square root 
of the expectation value of p' [13] , 

3) If both of states are diagonal in the same basis, such 
as the thermal equilibrium state, the fidelity (or classical 
fidelity) can be calculated as 



F{p,p') 



E 



(9) 



Example: 

If a spin is coupled to environment, it can be 
described by a reduced-density matrix. For 
two reduced-density matrices 



a 
b 



P 



c 
d 



the fidelity can be calculated as 
F{p,p') = y^+Vbd. 



Though the fidelity itself is not a metric, it can be used 
to define a metric on the set of quantum state, i.e. 

0b{p,p') = co s-' [F{p,p% (10) 

Db{p,p') = V2-2f(p,p0, (11) 

Ssip.p') - yWFWOP, (12) 

and 



called commonly as Bures angle, Bures distance [2£ 
sine distance [30j, respectively. 

Besides the above well-accepted definitions, there are 
some alternative definitions of the fidelity. For example, 
Chen et al [si'l proposed 

|2 



\F{p,p') 
1-r 1 



(13) 



tripp') + Vl-tr(p2)^l-tr(p'2) 



2 2 

where r = l/{d — 1) with d being the dimension of the 
system. This definition has a hyperbolic geometric inter- 
pretation, and is reduced to Eq. [8] in the special case of 
d = 2. The definition [Eq. (fT3|) ] was recently simplified 
to [33] 



]F(p,p')]' = tr(pp') + ^1- tr(p2)^l _ tr(p'2). (14) 

Obviously, one of advantages of the above definitions is 
that the fidelity can be easily evaluated for arbitrary 
mixed states. Nevertheless, it seems that for two density 
matrices of two sets of mutually independent events, Eq. 
(|14p gives a nonzero value. Therefore, another definition 
of the fidelity was proposed [13], i-e. 



F{p,p') 



Mpp')\ 



Vtr(p2)tr(p'2) 



(15) 



The fidelity has been widely used in many fields. In 
quantum information science, the fidelity between quan- 
tum states have been proved useful resources in ap- 
proaching a number of fu ndamental prob lems such as 
quantifying entanglement [llll Ill2l . Ill3l ]. There are 
also many interesting works on the fidelity in adia- 
batic processes. For example, the adiabatic fi delity was 
used to describe atom-to-molecule conversion |ll4 IllSj ] 
in atomic systems and the time evolution in a Bose- 
Einstein condensat eillGl . IllTl . IllSl . Physicists working 
on quantum chaos Il2l| use quantum fidelity 

(Loschmidt echo) to measures the hypersensitivity to 
small perturbations of quantum dynamics. In the latter 
case, the fidelity usually depends on the time. Interested 
readers can find more details about the Loschmidt echo 
in a recent review article by Gorin et al 32]. In the fi- 
delity approach to quantum phase transitions, which will 
be introduced in this review, the fidelity depends on the 
adiabatic parameter (or driving parameter) of the Hamil- 
tonian, and is usually static. 



B. Quantum state overlap and adiabatic evolution 

To well understand the fidelity in the ground-state 
state evolution, in this subsection, we take a 1/2 spin sub- 
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FIG. 1: (Color online) A single 1/2 spin state defined on a 
Bloch sphere. 



There are two parts in the overlap. The real part de- 
notes the difference in the geometrical structure, while 
the imaginary one corresponds to the overall phase dif- 
ference. 

Though the overlap shows also the similarity between 
two states, it does not show the response of the state 
at a given point to the adiabatic parameter (p. For this 
purpose, we expand the overlap around a given </> as 



a2 



where 



*(6 



2' 



(21) 

(22) 
(23) 



The linear term is the Berry adiabatic connection, which 
contribute a Pancharatnam-Berry phase [l22L Il23l | to the 
spin as the magnetic field rotates adiabatically around 
cone direction (the dotted circle in Fig. [1]), i.e. 



jected an external magnetic field as a warm-up example. 
The Hamiltonian of a free spin under an arbitrary field 
B(i3 sin 6 cos 0, B sin 9 sin 0, B cos 0) with magnitude B is 



H 



B cr, 



(16) 



where a = [a^ ,a'^ , cr^) are the Pauli matrices. In basis 
I t), I i), Pauli matrices take the form 



cr = 



1 

1 



, cr 



y 



-i 

1 



, cr 



1 

-1 



(17) 



Then the Hamiltonian (I16p matrix can be rewritten as 
cos 9 e~"^ sin 



H 



-B 



sm ( 



cost 



(18) 



The Hamiltonian can be easily diagonalized and the 
spin's ground state, with eigenenergy Eq = — i?, is 

|vI/(0,0))=cos^|r)+e'^sin^U). (19) 

Here 9 and (j) define a point on the unit three-dimensional 
Bloch sphere to which the spin points to(see Fig. [1]). The 
state can be multiplied by an arbitrary global phase. Ob- 
viously, 9 and (p can be regarded as adiabatic parameters. 
For simplicity and without loss of generality, we fixed 9 
first. Then the overlap between two states corresponding 
to two points on the ring of a given 9 is 



f{9,<P;9,<p') = {^{9,c^)\^{9,<f,')) 



cos 2 + ^ 



(20) 



7(0,0) 



d 



,sin^-. 



(24) 



The phase equals to the solid angle of the cone if the spin 
rotates one periodicity. The Berry connection must be a 
purely imagnary number because of 



= 0. 



(25) 

The global phase can be rectified by a gauge trans- 
formation e~'^'^^^''^\ which can compensate the geometric 
phase "f{9, </>) accumulated during the adiabatic evolu- 
tion. The new state becomes 



cos-|T) + e^*sin-U; 



(26) 



Then the overlap between two geometrically similar 
states becomes 



fi9,(P;9,(P + S<P) = exp 



X I cos^ - 4- e-'^'l' sin^ ■ 



W 2 .26* ^9 
' %\VL - cos"^ - 
2 2 2 



•(.27) 
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The most relevant term is then the second derivative of 
the overlap. Moreover, the gauge transformation not only 
eliminates the Berry adiabatic connection, but also mod- 
ifies the second order term. After the phase rectification, 
the second-order term is reduced to a minimum. On the 
other hand, the phase rectification denotes mathemati- 
cally a rotation in the complex plane, which makes the 
overlap be a purely geometric quantity. 

If we take the modulus of the overlap, it becomes the 
fidelity 



Iff 



>sin2- 



(<5</>)2 . 2 e 

— — sm - 

2 2 



f-(<50)2sin2-cos2- 



(28) 



Then the fidelity, if we express it in a series form, becomes 
F=|/| = l-M!3in^^cosi-H.... (29) 



Therefore, the leading response of the fidelity to the adi- 
abatic parameter is its second derivative. This is quite 
natural because the fidelity can not be large than its up- 
per limit 1, it must be an even function of the pertur- 
bation of the adiabatic parameter. The leading term is 
called fidelity susceptibility in some literatures because 
it is physically a kind of structure of the driving term. 





FIG. 2: (Color online) A sketch of a ground-state level cross- 
ing (LEFT) and the first-excited state level-crossing(RIGHT) 
as the system's driving parameter varies. 




FIG. 3: (Color online) A sketch of a quantum phase transition 
occurred at Ac (LEFT) and corresponding expected behavior 
of the fidelity F = {*o(Ai)l*o(A2)) as a function of A = 
(Ai + \2)/2 for a fixed 5A = A2 - Ai (RIGHT). 



XF 



siv? — cos^ - 
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(30) 



Though the phase rectification can change the Berry adi- 
abatic connection and the second derivative of the over- 
lap, the fidelity susceptibility does not change. This phe- 
nomenon is due to the simple reason a gauge transforma- 
tion cannot affect the modulos of the overlap. 

On the other hand, when we study quantum phase 
transitions occurring in a quantum-many body system, 
the ground-state wavefunction is usually defined in the 
real space, then the imaginary part of the overlap does 
not appear. If the adiabatic parameter is defined on the 
flat manifold, the linear correction is zero. The second 
term is the most important, ft denotes the leading re- 
sponse of the wave function to the adiabatic parameter. 
Though for the present case it is simply a constant due to 
the rotational symmetry of (f), it might become singular 
for a many-body system in the thermodynamic limit. 

Now we consider another case of fixing both 9 and 
(j), and changing the magnitude of the external field. If 
B > 0, the ground state is 



with eigenenergy B. A ground-state level-crossing oc- 
curs at the point B = 0. Then the fidelity shows a very 
sharp drop a,t B — due to the level-crossing between 
two orthogonal states. While if we expand the fidelity in 
term of B, one may find that either the fidelity suscep- 
tibility or the Berry adiabatic connection is zero except 
for B = 0. The point is a singular for both of the fi- 
delity susceptibility and the Berry adiabatic connection. 
In many studies on the Pancharatnam-Berry phase, this 
level-crossing point is regarded as a monopole in the pa- 
rameter space. 



III. FIDELITY AND QUANTUM PHASE 
TRANSITIONS 



The fidelity and its leading term introduced in the last 
section is illustrative. In this section, we try to establish 
a bridge between quantum phase transitions and the fi- 
delity in considerable detail through the one-dimensional 
transverse-field Ising model. 



|*(0, </>))= cos- 1 r)+e^^ sin -li), (31) 

with eigenenergy —B, while if _B < 0, the ground state 
becomes 



|v|/(0,</.))=e^^sin-|T)-cos-U), (32) 



A. Quantum phase transitions: fidelity perspective 

Without loss of generality, the Hamiltonian of a gen- 
eral quantum many-body system, which might undergo 
a quantum phase transition in parameter space, can be 
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written as 

H{\) = Ho + XHi, (33) 

where Hj is the driving Hamiltonian and A denotes its 
strength. According to quantum mechanics, the system 
satisfies the Schrodinger equation 

iJ(A)|*„(A)> = £;„|*„(A)), n = 0, 1, • • • , (34) 

where En is the eigenenergy and set to an increasing 
order Eq < Ei < E2 ■ ■ ■ , and |vI'„(A)) defines a set of 
orthogonal complete bases in the Hilbert space, i.e. 

^|vI/„(A))(vl/„(A)|=/. (35) 

n 

As the driving parameter A varies, the energy spec- 
tra are changed correspondingly. The quantum phase 
transition occurs as the ground-state energy undergoes a 
significant change at a certain point. Precisely, its first- 
or higher-order derivative with respect to the driving pa- 
rameter becomes discontinuous at the transition point. 
There are two distinct ways. The first one is the energy 
level-crossing occurring in the ground state (left plot of 
Fig. [2]). The second is that the level-crossing occurs 
only in the low- lying excitations 0] , and the ground state 
keeps nondegenerate (right plot of Fig. (2). For both 
cases, the structures of the ground-state wavefunction 
become qualitatively different across the transition point. 
That is, if we compare two ground states on both sides 
of the transition point, their distance is very large; while 
if we compare two ground states in the same phase, their 
distance is relatively small. Therefore, if we calculate 
the fidelity between two ground states, i.e., the fidelity 
of 5*0 (Ai) and ^'o(A2) at two slightly separated points 
Ai(2) with fixed S = Xi — A2, it should manifest a mini- 
mum at the transition point, as shown in Fig. [3l Such 
a fascinating perspective for quantum phase transitions 
was first observed in the one-dimensional transverse-field 
Ising model [35ll3^. 

Obviously, the fidelity between two ground states does 
not bear any apparent information about the difference 
in order properties between two phases. Instead, it is a 
pure geometric quantity of quantum states. In its ap- 
proach to quantum phase transitions, one of obvious ad- 
vantage is that no priori knowledge of order parameter 
and symmetry-breaking is required. For example, if a 
quantum phase transition is induced by the ground-state 
level crossing, then the two crossing states at the transi- 
tion point are orthogonal, then the overlap between them 
is zero; while the fidelity almost equals to one in other 
region away from the crossing point. Therefore, it is be- 
lieved that the fidelity can describe quantum phase tran- 
sitions in its own way. 

B. Example: the one-dimensional transverse-field 
Ising model 

T he on e-dimensional transverse-field Ising model [l25l 

Il26l . Il27| is one of the simplest models which can be 




FIG. 4: (Color online) A sketch of the one-dimensional 
transverse-field Ising model with periodic boundary condi- 
tions. Two arbitrary neighboring spins interact with each 
other by the Ising interaction ajaj^i (dotted fine). All spins 
are subject to an external field h along z direction. 
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FIG. 5: (Color online) A schematic ground-state phase dia- 
gram of the one-dimensional transverse-field Ising model. At 
the left hand side of he, the ground state is in a ferromagnetic 
long-range order phase whose low-lying excitations are flipped 
spin quasi-particle. While at the right hand side of he, the 
ground state is a fully polarized phase whose low-lying exci- 
tations are domain wall quasi-particles. 



solved exactly [l28l . Il29l | in the field of condensed matter 
physics. Due to its simplicity and clear physical pictures, 
the model is often used as a starting model to test new 
physical ideas and approaches, among which the fidelity 
does not make an exception. The following procedure is 
standard, and the final expression of fidelity is obtained 
by Zanardi and Paunkovic [36| . 

The Hamiltonian of the one-dimensional transverse- 
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field Ising model with periodic boundary conditions reads 



H 



N 

X 



(36) 
(37) 



where h is the transverse field and N is the num- 
ber of spins. As inferred from the model's name, the 
Hamiltonian describes a chain of spins with the nearest- 
neighboring Ising interaction along x-direction, and all 
spins are subject to a transverse magnetic field h along 
the z-direction (Fig. |4]). 

The Hamiltonian is invariant under translational op- 
eration. Moreover, unlike usual spin systems, the z- 
component of total spins in this model, i.e. 



N 



(38) 



is not conserved. Instead, if we introduce 

a+ = i (a- + lay) , a- = i (a- - za^) , (39) 



one-dimensional Ising model. Defining the eigenstates of 
cr^ as 

l-> = ^(IT) + U)),IO = ^(IT)-U)), (45) 

the doubly degenerate ground states of the Hamiltonian 
take the form 



l*i> 

l*2> 



(46) 
(47) 



which are of ferromagnetic order. The ground-state prop- 
erties change as the external field h turns on. Because 
of 



(48) 



the magnetic field mixes l^*!) and \^2) and the ground 
state becomes non-degenerate for a finite system. De- 
spite of this, the ground-state property does not change 
qualitatively. The ground state still manifests the fer- 
romagnetic long-range order. Precisely, the correlation 
function 



and 















a = 


(Jo) 




T) = 




U) = 1 T 



(40) 



(41) 



the Hamiltonian p6p can be transformed into 

N 



(44+1 + '^j+O 



(42) 



then we can see that the off-diagonal terms in the Hamil- 
tonian (|42)) either exchange the state of a pair of anti- 
parallel spins, or flip two upward spins to downward or 
vice versa. So they do not change the parity of the sys- 
tem. This property defines a classification of subspaces 
based on the parity operators, i.e. 



N 



(43) 



and the Hamiltonian cannot change the parity of the 
state, i.e. 



[H, P] = 0. 



(44) 



Therefore, we have two subspaces corresponding to parity 
P = ±1 respectively. 

The ground state of the one-dimensional transverse- 
field Ising model can be understood from its two limiting 
cases. If /i = 0, the Hamiltonian becomes the classical 



(vI/ok,-aJ+,,|vI/o) - (*okJ|*o)(*okJ+,.|*o>, (49) 

does not vanish even lir oo. The correlation function, 
therefore, can be used as an order parameter to describe 
the phase in the small h region. While if /i — > oo, the Ising 
interaction is neglectable, all spins are fully polarized 
along z-direction. The ground state is non-degenerate 
and takes the form 



l^-o) = I TTT •••??)• 



(50) 



In this limit, the correlation function Eq. does 
not show long-range behavior. Therefore, a quantum 
phase transition between an ordered phase to a disor- 
dered phase is expected to occur as h changes from zero 
to infinite. A schematic ground-state phase diagram of 
the model is shown in Fig. [Sj 

In order to discuss the fidelity in the ground state, we 
now diagonalize the Hamiltonian in detail. We need three 
transformations, i.e., the Jordan- Wigner transformation 
(l30j |. Fourier transformation, and Bogoliubov transfor- 
mation. 

The Jordan- Wigner transformation: The Jordan- 
Wigner transformation maps 1/2 spins to spinless 
fermions, that is 



exp 



exp 



n-l 

in c]cj 

n-l 
'ITT > CjCj 



n 



cr|c„, (51) 



n 



^]cl (52) 



1 2c]jCn, 



(53) 
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where cj^ and c„ are fermionic operators and satisfy the 
anticommutation relations 



(54) 
(55) 



After the Jordan- Wigner transformation, the Hamilto- 
nian becomes 



JV-l 



H = - 



,t J 



{c\cN + c^fCi j exp 



cicjv 1 exp 



AT 



N 



The exponential factor 
P = exp 



(56) 



N 



3 = 1 



is nothing but the parity of the system which is a con- 
stant, i.e. for periodic boundary conditions P = —1 and 
antiperiodic boundary conditions P = 1. The Hamilto- 
nian can be simplified as 



H = - 



N 

mi 



Cj + l + c 



i=i 

N 



(44+1 + Cj + lCj) 



1 - 24c, 



(57) 



Fourier transformation: Since the Hamiltonian is in- 
variant under translational operation, we can perform 
standard Fourier transformation. For the present case, 
the transformations are 



~ikj 



(58) 



where the momentum fcs are chosen under conditions: 

k=< ^_ \ , (59) 



2n7r 

N 



with n = 0, 1, 2, • • • — 1. Then the Hamiltonian can be 
transformed into fc-space form, 



- ^ 1^(2 cos k - 2h)clck + i sin k (c^ j,c| 

k 

-Nh. 



f. + C^kCk 

(60) 



Bogoliubov transformation: Obviously, the quadratic 
Hamiltonian can be further diagonalized under the fa- 
mous Bogoliubov transformation: 



Ck 



C-k 

elk 



Ukbk + ivkbl^., 
Ukb\ - ivkb-k, 
Ukb-k - ivkbl, 
Ukbik + ivkbk, 



(61) 



where bk and fcj, are also fermionic operator and satisfy 
the same anticommutation relation as Ck and cj, . Because 
of this, one can find the coefhcients in the transformation 
([61]) should satisfy the following condition 



Vk 



-v-k, Uk+'"k 



1. 



So we can introduce trigonal relation 

Vk = sin6'fc, Uk = cos 6k. 



Inserting the Bogoliubov transformation into Eq. (|60|). 
the coefficients are determined by 



(62) 



(63) 



cos 29k — 
sin 20fc = 



cos(/c) — h 



v/l-2/icos(fc)-)-/i2' 

sin(fc) 
^l-2/icos(fc) + /i2^ 



(64) 



such that the Hamiltonian becomes a quasi-free fermion 
system. 



where 



e(fc) yj\ - 2hcos{k)+h'^ 



(65) 



(66) 



is the dispersion relation of the quasi particles. The dis- 
persion relation shows that the thermodynamic system 
is gapless only at h = I, and gapped in both phases of 
< h < 1 and h > 1. Therefore, the quantum phase 
transition occurs at the point h = 1. 

The ground state: The ground state of the model is 
defined as the vacuum state of bd'^o) = where 



bk = cosOkCk - ism OkC_i., 
b-k = cos0fcC_/c -I- jsin^fcc^. 

Since the condition 

a|0)fc|0)_fc + 6|l)fc|0)_fe 



gives 



' +c\0)k\l)-k + d\l)k\l)-k 



a = cos Ok, b = 0,c— 0,d = i sinOk, 



(67) 



(68) 



(69) 
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FIG. 6: (Color online) The fidelity of tfie one-dimensional 
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FIG. 7: (Color online) The fidelity of the one-dimensional 
transverse-field Ising model for various Sh. Here A'^ = 290. 



the ground state takes the form 

l*o(M) = n (cos0fc|O)fc|O)_fc-Hisin0fc|l)fe|l)_fe). 

fe>0 

(70) 

The low-lying excitation can be obtained by applying fcj, 
to the ground state. In the ferromagnetic phase, the 
excitation is visualized as a quasi-particle of a flipped 
spin; and as a domain wall quasi-particle in the fully 
polarized phase (Fig. [5]). 

Fidelity: Once the ground state is obtained expl icitly , 
the fidelity between h and h' can be calculated as |l3Cll | 



because the ground-state wavefunction (iTOl) is already a 
product state. 

Fig. [Hlshows the ground-state fidelity of the transverse- 
field Ising model as a function of h with parameter dif- 
ference Sh = 0.01. The numerical results of a smaller 
sample, say 20 sites, have also been compared with exact 
numerical computations and the agreement is essentially 
perfect (see Table II of section IVIip . As expected, the 
quantum critical region is clearly marked by a sudden 
drop of the value of fidelity. The behavior can be ascribed 
to a dramatic change in the structure of the ground state 
of the system during the quantum phase transition. The 
drop becomes sharper and sharper as the system size in- 
creases. Meanwhile the fidelity in the non-critical region 
is also reduced though the reduced magnitude is smaller 
than that at the critical point. This property can be in- 
terpreted due to the increasing of number of degree of 
freedom. Actually, in the thermodynamic limit, the fi- 
delity between two different ground states might be zero, 
no matter how small the difference in parameter Sh is. 
That is the two ground states are orthogonal to each 
other. This phenomena has been studied in quantum 
many-body systems, a nd is known as the Anderson or- 
thogonality catastroi)he \l2'^ . Fig. [7] shows the fidelity 
for a given size system but various Sh. The figure is easy 
to be understood. The larger the distance between two 
points in the parameter space, the larger the distance 
between the two corresponding ground states. 



IV. FIDELITY SUSCEPTIBILITY, SCALING 
AND UNIVERSALITY CLASS 



The leading term of the fidelity and dynamic 
structure factor 



The differential form: A sudden drop of the fidelity 
caused by the ground-state level-crossing is too obvious 
to be interesting enough. People are interested in those 
ground-state wavefunctions which are differentiable in 
parameter space. Therefore, the overlap between two 
ground states at A and X + 6X can be defined as 



/(A,A + 5A) = (^o(A)|*o(A + <5A)). 



(72) 



F{h,h) — \{'i>o{h)\'i>o{h))\ — Y\_'^os{Ok Ok)- (71) Performing series expansion, the overlap becomes 



k>0 

As is always emphasized, the fidelity is purely a geomet- 
rical quantity since it is an inner product between two 
vectors. Eq. (|7ip refresh our mind on this point because 
the expression is just an angle between two vectors. It is 
also consistent with the fourth fidelity axiom of Eq. ^ 



f{\,X + SX) = l + 5A^*o(A) 



*o(A) 



d 



(73) 



— vI/o(A)\-f..(74) 
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The fidelity, as the absolute value of the overlap, then 
becomes 



the p arenthesis of Eq. ((8T|l is called quantum metric 
tensor jlOOl Ill0l | or the Ricmann metric tensor (4J] 



\fiX,X + SX)f 



= 1 + SX 



<SX) 

1 



*0 

d 



dX 
9A2 



d_ 
dX 



dX 



*o ) ( *o 
1 / 

dX^ 



dX 



•(75) 



The linear correction must be zero. There are two rea- 
sons. The first is due to the normalization condition, 
i.e. 



d_ 

dX 



d_ 
dX 



*o ) = ^(^-ol^o) = 0. (76) 



The second is that the |/| must be small than 1, then the 
leading term must be an even function of 5X. Therefore, 



F{X,X + 5X) = \-^-^XF 



(77) 



where xf denotes the fidelity susceptibility of the ground 
state. 



1 / d 

2 \dXa 
d 

dXa 



d 
dXb 



1 



d 



2 \ dXb 
d 



d 
dXa 



dXb 



(82) 



The quantum metric tensor is symmetric under exchange 
of the index a and b. It is the re al pa rt of a more general- 
ized quantum geometric tensor [l09j | of the ground state. 
Precisely, if we defined the projection operator 



P 



(83) 



which projects out the ground state, the quantum geo- 
metric tensor then is defined as 



P 



d 

dXb 



Therefore, 



gab — Re Tab, 



(84) 



(85) 



XF 



d_ 
dX 



d_ 

dX 



d ^ 

dx ° 



9 T 

dx ° 



(78) 



This is the differential form the fidelity susceptibility. 

On the other hand, if the ground-state wavefunction 
is defined in the multi-dimensional parameter space, say 
A — {Xa], a = 1, 2, • • • A, the overlap between two states 
at A and A' = A (5A is 



f{X,X + 5X) = l + Y,5Xal-^o 



_d_ 

dXa 



ah 



d d 



dXa dX\ 



(79) 



The fidelity susceptibility becomes 
^ dx dX I 2 \dXa 



dX, 



1 / d 

2 \9A^ 
d 

dXa 



d 



dXa 



9A, 



(80) 
(81) 



and the imagnary part of Tab defines a 2-form phase. 



Vab = 2Imraf, 



(86) 



The 2-form phase Vab plays a very important role in ge- 
ometric phase. Its flux gives the Berry phase. While the 
quantum geometric tensor provides a natural means of 
measuring distance along the evolution path in parame- 
ter space. The distance between two ground states can 
be expressed in the differential-geometrical form, i.e 



ds^ = y^^gabSXgSXb. 



(87) 



In addition, if the ground state of the system evolves 
adiabatically from A to A' at a given path S, the quantum 
distance Rq in the parameter space is 



Rq= f y^gabdXadXb. 
^^^V ab 



Therefore, if we do geodesies, 



dRq = 0. 



(89) 



where the vector dXa/dX denotes the direction of the 
short displacement SX in parameter space. The term in 



we can in principle find the shortest path connecting the 
two ground states at A and A'. 
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Example: 

Take the spin in an external field as an exam- 
ple, its ground state is 

|*(^,<^))=cos^|T)+e'*sin^U). 

The quantum metric tensor takes the form 
1 



g 



sm( 



ds^ = de'^ + sin edd)'^ 



so 



which is the metric on the sphere of parame- 
ters. 



The perturbation form: We concern mainly on the 
fidelity in continuous phase transitions. That is, the 
ground state of the Hamiltonian is nondegenerate for a 
finite system. Therefore, as the point A -I- (5A closing to 
A, the ground-state wavefunction can be obtained, to the 
first order, as 



|*o(A + (5A)) 
where 



(90) 



(91) 



is the hoping matrix of the driving Hamiltonian Hj. 
Therefore, if we normalized the wavefunction 1^*0 (A -I- 
SX)), the fidelity becomes, to the leading order. 



1 - 5X 



| (vI/„(A)|g,|vl/o(A))| 

[En{X)-EoiXW 



(92) 



Obviously, the fidelity depends both on A and SX. The 
most relevant term is the leading term in Eq. i.e. 
the second order derivative of the fidelity with respect to 
SX. The term actually defines the response of the fidelity 
to a small change in A. The fidelity susceptibility can be 
obtained as 

Xf{X) = hm 



<5A^0 (5A2 



disxy 



With Eq. dUD, it can be rewritten as [isl.lii 



Xf(A) 



E 

n#0 



\{^„{X)\Hj\MX))\' 
K(A)-i?o(A)]2 



(93) 
(94) 



(95) 



This is the summation form of the fidelity susceptibil- 
ity. The form establishes a relation between the struc- 
ture difference of two wavefunctions and low-lying energy 
spectra. 



Example: 

To understand Eq. (pSj). we still take the spin 
in an external field [Eq. (|16p ] as an example. 
The driving term in the Hamiltonian at a fix 
point can be obtained as 



dH 



iB 



e-"P sin ( 
-e*"^ sin 61 



The excited state is 

|*i(0,0))=e'^^sin^|t)-cos^U) 

Then the hoping matrix between the ground 
state and excited state takes the form 



dH 



= iBe'^ sin( 



The energy different between two state is 2B, 
then the fidelity susceptibility becomes 



1 



Xi^ = - sm 



which is the same as Eq. (pO|) . 



On the other hand, according the perturbation theory, 
the second order perturbation to the ground-state energy 
takes the form 



|(^„(A)|g,|vl/o(A))| 
Eo{X)-E„{X) 



(96) 



Obviously, Eq. (|95|) and Eq. (I96|) are very similar in 
their form except for different exponents in both denom- 
inators [HH. Therefore, one might expect that the origin 

(2) 

of the singularity of the fidelity susceptibility and i?g 
are both due to the vanishing of the energy gap though 
the fidelity susceptibility shows a sharper peak than Eq ' . 

(2) 

For the finite-order phase transition, however, Eq can 
be still a continuous function of the driving parameter, 
then there is no reason to require that the fidelity suscep- 
tibility shows singular behavior in high-order(> 2) quan- 
tum phase transitions. It was also pointed out later that 
the fidelity susceptibility might be related to the third 
energy perturbation ^63,] , 



1 



E, 



(3) 



XF 



^ [E, 



Eom-Eo) £;(!)■ 



(97) 



Therefore, the fidelity susceptibility might not be able to 
witness those phase transitions of infinite order [i^ [Blj , 
such as the Beresinskii-Kostcrlitz-Thouless transition. 

The fidelity susceptibility as a kind of fluctuation: The 
hoping matrix (\E'„|i//|^'o) implies dynamics behaviors of 



14 



the fidelity susceptibility. Similar to the linear response 
theory, one can define the dynamic fidelity susceptibility 



XF 



M = E 



\{^n\Hj\^0)\' 

[En - Eo]^ + uj^ ■ 



(98) 



Performing a Fourier transformation, the dynamic fi- 
delity susceptibility becomes 



^ En — Eq 



(99) 



The energy difi^erence in the denominators can be can- 
celed if one take a derivative with respect to r, the dy- 
namic fidelity susceptibility then is 



-ttGi{t)9{t) +TrGii-T)ei-T). 




(100) 



G/(r) = (*o|i?/(T)i?/(0)|«'o)-(«'o|i?/|*o)'(101) 
Hi{t) = e^(^)^i^/e-■f^(^)^ (102) 

with T being the imaginary time. Performing an inverse 
Fourier transformation, we can obtain 



Xf(w) 



1 



Sm(LUT)Gl{T)dT. 



The fidelity susceptibility then becomes 
1 /"^ 

Xf = lim — / sm{uJT)G I {T)dT. 
u — >o uj Jq 



(103) 



(104) 



For any finite system, the correlation function Gi{t) de- 
cays in the large r limit, the above limit satisfies the 
Lebesgue's convergent theorem, the fidelity susceptibil- 
ity, finally, has the form [isj 



/•oo 

XF = TGl{T)dT. 

Jo 



(105) 



Therefore, the fidelity susceptibility is nothing but a kind 
of dynamics structure factor of the driving Hamiltonian. 

The Eq. (|105p is remarkable because it connects the 
fidelity to dynamical response of the system by the driv- 
ing Hamiltonian Hj. In this way, the adiabatic evolution 
of the ground state and the fidelity susceptibility are ex- 
pressed in terms of standard quantities in linear response 



theory and their physics content is clarified. Tradition- 
ally, quantum phase transitions are said to be driven by 
quantum fluctuations, which originate from the Heisen- 
berg uncertainty relation. In the Hamiltonian 



(106) 



The the second order perturbation to the ground-state 
energy 



E, 



(2) 



[(*o|ff/(r)iJ/(0)|*o) - (*o|ff/|*o)'] dr 



(107) 



is also a kind of fluctuation. Clearly, both Eq and xf 
become zero if [Ho,Hi] = 0. Therefore, the expression 
(|105p provides a new angle to understand the role of 
quantum fluctuation in quantum phase transitions. 

The quantum metric tensor: In case that the Hamil- 
tonian is defined in a high-dimensional parameter space, 
the fidelity susceptibility becomes a metric tensor. The 
Hamiltonian in A-dimensional parameter space reads. 



H = yXaHa, a = l,...,A, 



where AqS are coupling parameters. Clearly 

dH 



dXa' 



(108) 



(109) 



In the parameter space, we can always let the ground 
state of the system evolves along a certain path, i.e 



— Aa(A), 



(110) 



where A plays a kind of driving parameter along the evo- 
lution line. Therefore, the driving term in the Hamilto- 
nian at a given point A is 



dH ^dX 



dX 



(111) 



Then the fidelity susceptibility along this line can be cal- 
culated as 



Xf(A) 



|(*„(A)|H,|vI/o(A))|^ 



[En{X)-E,{X)Y 
^ dXa dXb 

ab 



TrOn TTnO 

Z_. dX dX Z^jEn{X)-Eo{X)]^''- ' 



where if™,") = (*„,(A)|i?„(b) |*„(A)). So the quantum 
metric tensor takes the form 



9ab = 



TjOn TjnO 



[En{X) ~ Eo{X)Y 



(113) 



This is the perturbative form for the quantum metric 
tensor. 



Oab 



TGab{T)dT, 



(114) 



15 



where Gab{T) is a time-dependent correlation function, 

Gabir) = 0(T)((*o|i^a(T)i/b(O)|*o) 

-(*o|i^a(0)|*0>(*o|i^6(0)|*o))- (115) 

Clearly, Eq. (|114p is an extension of Eq. (jlOSp . Both of 
them clarify the physics content of the fidelity approach 
and the role of quantum fluctuation in the adiabatic evo- 
lution. 



B. Scaling analysis and universality class 

In physics, if a physical quantity depends linearly on 
the system size, the quantity is extensive (additive) , such 
as energy; and if it is independent of the system size, it 
is intensive, such as energy density. If we define A as 
the energy gap between the ground state and the lowest 
excitation with nonzero (^„(A)|-ff/|^'o(A)), the fidelity 
susceptibility satisfies the following inequalities 

1 ^\{^^i\)\Hj\^o{\))? 1 d'EiX) 



K(A)-£;o(A)] 2 A aA2 ' 



< 



^Y.\{m^{\)\Hi\^,{\))\ 



(116) 



-^[(vl/o(A)|i/,2|*o(A)) - {^,{X)\H:\^,{\))^] 



The inequalities tell us some useful information about 
the properties of the fidelity susceptibility away from the 
critical point: 

1. If the system is gapped and Hi behaves like a single 
particle, xf is an intensive quantity. 

2. If the system is gapped, the fidelity susceptibil- 
ity share the same (or smaller) dependence on the 
system size as the second order derivative of the 
ground state energy. 

3. If Xf is a superextensive quantity, the ground state 
of the system should be gapless. 

Therefore, unlike conventional physical quantities, the 
fidelity susceptibility can either be an extensive quan- 
tity, like the energy, or has other type of dependence on 
the system size. To have an explicit view of the depen- 
dence of the fidelity susceptibility on the system size, it is 
useful to perform scaling transformation, which is firstly 
carried out by Venuti and Zanardi [i^l- Without loss of 
generality, the interaction Hamiltonian can be written as 
a summation of local terms, i.e. 



and number of site 



N 



(117) 



(118) 



where d is the real dimension of the system. Then the 
fidelity susceptibility becomes [i^ 



g=5:|rC(r,r)dr, 



(119) 



where 



C(r,T) = (1'o(A)|y(r,T)y(0,0)|^'o(A)) 
-(vI/„(A)|F(r,0)|vI/„(A))(vI/„(A)|y(0,0)|vI/o(A)). (120) 

In the vicinit y of the critical point, the scaling transfor- 
mation goes |l3ll Il32l | 



s^T, V{r') = s-'^^'F(r), 



(121) 



where s > 1, C is the dynamic exponent, and Ay is the 
scaling dimension of V{r)^ one can find that 



1 



{L'/sY 
Therefore 



J T'C{r',T)dT'. 



x'f 



Xf 



{L'Y s<i+2C-2Av. Id ' 



(122) 



(123) 



which defines the scaling transformation of the fidelity 
susceptibility. Around the critical point, the correlation 
length is divergent and the only length scale is the system 
size itself, the fidelity susceptibility scales like 



Xf rd+2C-2Av 



(124) 



The expression is interesting. It establishes a relation 
between the size dependence of the fidelity susceptibility, 
the dynamic exponent, and the scaling dimension of the 
driving term. 

In most cases, we can regard Xf/L''' as an intensive 
quantity in general quantum phases because the system 
is usually gapped. Therefore, the singular behavior of 
the fidelity susceptibility is due to 



d + 2C > 2A 



V 



(125) 



at the critical point. If the fidelity susceptibility diverges 
as 



XL 



1 



|A-Ae|" 

the scaling exponents satisfy 

rf + 2C-2Ay 



(126) 



(127) 



where v is the critical exponent of the correlation l eng th. 
A similar relation was also obtained by Gu et al [49[ in 
their studies on the one-dimensional asymmetric Hub- 
bard model. 
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However, Eqs. (|124m27t are not universally true in all 
quantum phase transitions. Considering a general corre- 
lation function 



C(r,T) 



1 



„2A 



(128) 



the fidelity susceptibility, in the thermodynamic limit, 
becomes 



XF 



E 



-/(rri/«)dT, 



,,2Av 



cx 



V-2C 

InL. 



2Av -2C^d 



(129) 
(130) 
(131) 



Such an explict expression implies that the fidelity sus- 
ceptibility is not alway extensive. Even in a gapped 
phase, as we can see from the Lipkin-Meshkov-Glick 
model, it can be intensive. Therefore, if we define the 
size exponent above (below) the critical point as c?"*"(c?~), 
that is the fidelity susceptibility 



(132) 



Hence xf / L'^ is an intensive quantity in corresponding 
phases. Physically, since the fidelity susceptibility de- 
notes the response of the ground state to the adiabatic 
parameter and c?J is the dimensional dependence of the 
adiabatic evolution, has been called adiabatic dimen- 
sion. 

Therefore, the singular part of the fidelity suscepti- 
bility as an intensive quantity, around the critical point 
behaves Wgl 



XF 



1 



\X-\c\' 



(133) 



where denote the critical exponents of the fidelity 
susceptibility above or below the critical point. On the 
other hand, if the fidelity susceptibility around the criti- 
cal point shows a peak for a finite system, its maximum 
point at Amax scales like 



X(A = Amax) oc L'^ 



(134) 



where denotes the critical adiabatic dimension, and 
can be either analyzed from Eq. (|124p or obtained from 
numerical scaling analysis. The following function can 
include the above two asymptotic behaviors. 



X{\L) 



B{\ - K 



(135) 



where A is a constant, _B is a nonzero function of A around 
the critical point, and both of them are independent of 
the system size. According to Eq. (I135p . the rescaled FS 



is a universal function of the rescaled driving parameter 
L^{\ — Amax)j i-e.. 



XF(A)(A — Amax, L) - Xf(a)(A, L) 
Xf{\){\L) 



= /[L''(A-A„,ax)], 



(136) 

where v is the critical exponent of the correlation length. 
Then we have Ei, 



a 



(137) 



Therefore, unlike the second derivative of the ground- 
state energy, the fidelity susceptibility might have differ- 
ent critical exponents at both sides of the critical point. 
The above procedure is useful to determine the critical 
exponent from numerical computations, and has been 
used in some models. Nevertheless, it is still not com- 
plete. In some cases, the fidelity susceptibility shows 
logarithmic divergence around the critical point. Then 
we have a = , and Eqs. (|133m37|) should be changed 
correspondingly. For this case, at one side of the critical 
point, if xf / L'^'^ is intensive, and 



X(A = Amax) 



(X InL, 



(138) 



around the critical point. The logarithcmic divergence 
implies that 

1 - exp[Xi?(A) (A, L) - Xf(A) (A = Amax, L)] 

= /[L''(A-A,„ax)], (139) 

should be a universal function. 

Clearly, unlike conventional physical quantity that is 
either intensive or extensive, the fidelity susceptibility, 
as analyzed above, manifests distinct scaling behavior. 
This property might be due to both the relevance of 
the driving Hamiltonian under the renormalization group 
transformation and distinct dynamic exponent in differ- 
ent phases. 



C. Example A: the one-dimensional 
transverse-field Ising model 

The one-dimensional transverse-field Ising model gives 
us a simple example. For the Ising model, the fidelity is 

F{h,h') = \{^o{h')\^o{h))\ ^\{cos{eu~e'^). (140) 

fc>0 

At the point h, the fidelity susceptibility can be calcu- 
lated as 



XF^Y. 



k>Q 



where 



d9k 



sin k 



Mk ^ 1 

dX 2 1 + /i2 - 2/icosfc" 



(141) 



(142) 



17 



6 
5 

1 4 

<o 

3 3 

-a ^ 
E 

1 



_ ' 1 ' 


1 ' 


- 


- 

— >7V=70 




— A?=110 " 




--7V=150 - 




"iV= 190 







0.5 



1.5 




FIG. 8: (Color online) The fidelity susceptibility XfIN of 
the one- dimensional transverse-field Ising model for various 
system sizes. 



To find the scaling behavior of the fidelity susceptibihty, 
let us first consider the case of /i = 1, then 



Xf 



1 si^'^ ^ 

16 ^ (T 



(1 — cos A;)^ ' 



■n(N~\)lN 



tt/N 



fc2 



dfc, 



Therefore, we have 



(143) 

(144) 
(145) 

(146) 



for the one-dimensional transverse-field Ising model (see 
Fig. [8|). On the other hand, if /i ^ 1, the fidelity suscep- 
tibility becomes 



Xf 



N 



7r(JV-l)/JV 
/N 



sin k 



1 + /i2 - 2/i cos k 



dk. 



Obviously, there is no pole in the denominator of inte- 
grand. Xf then is an extensive quantity (see Fig. [8]), 



XF 

N 



sin(27ra;) 



1 + /i2 - 2/icos(27 



dx. 



(147) 



So the adiabatic dimenion for the Ising model is = 1. 
Explicitly, the integration can be evaluated by the residue 
theorem [tI]. Let 



sin(27ra;) = — 
2i 



cos(27ra;) = - z H — 
2 \ z 



Eq. (|147p becomes a contour integration along a unit 
circle. Then, we can obtain 



XF 

N 



1 



16(1 -/l2)' 



FIG. 9: (Color online) A sketch of the Lipkin-Meshkov-Glick 
model of 5 spins in which all 5 spins are mutually interact 
with each other (dotted lines) and subject to an external field 
h along z direction. 



for h < 1, and 



XF 

N 



1 



16/l2(/l2 - 1)' 

for /i > 1. At both sides of he — 1, we have 



= 1, 



(148) 



as the critical exponents. One can observe that in the 
both phases of the Ising model, the fidelity susceptibility 
is an extensive quantity. In addition, taking into account 
that the exponents correlation length — 1. The scaling 
relation = {d^ — d'^)/v is satisfied. 



D. Example B: the Lipkin-Meshkov-Glick model 

The one-dimensional transverse-field Ising model is not 
clear enough to explain the role played by the adiabatic 
dimension. For this purpose, we n ow take the Lipkin- 
Meshkov-Ghck model [Hi, [Hi, [Hil as an example. The 
Lipkin-Meshkov-Glick model was originall y intr oduced 
by Lipkin, Meshkov, and Click pM [iM [l35l | to de- 
scribe a collective motion in nuclei. The model consists 
of a cluster of mutually interacting spins in a transverse 
magnetic field A. Its Hamiltonian reads 



2XS, 



(149) 
7), (150) 



where S^. = J2i '^i f^i^ = Vi ^) spin 1/2 operators, 
S± = SxiiSy, and N the number of spins. The prefactor 
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FIG. 10: (Color online) The fidelity susceptibility in the 
ground state of the Lipkin-Meshkov-Glick model as a func- 
tion of A at 7 = 0.5. The inset denotes the scaling behavior 
of the maximum of the fidelity susceptibility, in which the 
sloDe of the line represents the size exoonent of the fidelity 
□ 




2 



FIG. 11: (Color online) The finite size scaling analysis is per- 
formed for the case of power-law divergence at 7 = 0.5 for 
system sizes = 2"(n = 12, 13, 14, 15, 16). The fidelity sus- 
ceptibility is considered as a function of system size and driv- 
ing parameter is a function of N^{\ — Amax) only, with the 
critical exponent u ~ 0.665 (From Ref. [13 )■ 



1 /N is to ensure finite energy per spin in the thermody- 
namic limit N 00. In the anisotropic case: 7 7^ 1, the 
Hamiltonian commutes vifith both the total spin S'^ and 
the parity P = Ylcrf, i.e. 

i 

[H,S']^[H,P]=0. (151) 

The symmetries significantly reduce the dimension of 
Hilbert space. For the present case, the ground state 
is a ferromagnetic when A > 0. Then the dimension of 
the subspace in which the ground state locates is N/2. 
This reduces the complexity of the problem, and one can 



study numerically a sample up to 2^^ spins using the 
standard diagonalization method. 

In the thermodynamic limit, the ground state of the 
system undergoes a second order quantum phase tran- 
sition at Ac = 1. If A > Ac, the system is fully mag- 
netized, while < A < Ac it is a symmetry-broken 
state. This conclusion was early drawn by the mean- 
field approaches JJ56s ,137.] . The finite-size scaling of this 
model was studied by the 1 /N expansio n in the Holstein- 
Primakoff single boson representatio n |138| and by the 
continuous unitary transformations |139l . [MOj. Recent 
studies also reveals a rich structure of the spectrum , fou r 
regions are distinghished in the parameter space [l4ll |. 
Besi des, peopl es also found that entanglement properties 
pi [l4a . [mI [list and fidelity [13 in the ground state 
of model can provide us a deep understanding on the 
quantum phase transition. 

The following results based on the fidelity approach 
done by Kwok et al [57]. Fig. (TU] show the depen- 
dence of the fidelity susceptibility on the driving param- 
eter for various system sizes. As expected, the peak 
around the critical point becomes sharper and sharper 
as the system size increases. Numerical analysis show 
that xf(A = Amax) ~ N< with dl ~ 1.33. The second 
observation is that the fidelity susceptibility shows differ- 
ent dependence on the system size in both phases. In the 
symmetry-breaking phase, X-f(A) ~ N, while in the clas- 
sical phase, xfW in an intensive quantity. Then adia- 
batic dimension of the fidelity susceptibility takes d~ = 1 
and = respectively. On the other hand, according to 
the scaling analysis discussed above, the rescaled fidelity 
susceptibility 

Xf(A^ Ainax,A^)"XF(A,A^) 

Xf{\,N) 

should be a function of N''{X — Xmax)- Fig- fTTl shows 
this function for various system size. All lines fall onto 
a single line for v ~ 0.665, Therefore, the exponents for 
Xf{X,N)/N'^ around the critical point, as an intensive 
quantity, 

± / 1/2 A < 1 

Therefore, the Lipkin-Meshkov-Glick model shows dis- 
tinct critical exponents around the critical point. Fortu- 
nately, the Lipkin-Meshkov-Glick model is also an exactly 
solvable model. The exact results then help us the check 
the numerical results obtained by scaling analysis. 

In the region A > 1, the Lipkin-Meshkov-Glick model 
can be dia gona lized using the Holstein-Primakoff trans- 
formation (l38| 

= S'-a^a, (153) 
S+ = {2S - a^a)^^^ a, (154) 

where a(a^) is bosonic annihilation(creation) operator 
satisfying [a, a^ = 1. In the large S{N) limit, the Hamil- 
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TABLE I: Critical exponents dj^, i^, d*, a* in various quantum phase transitions. The scahng relation a* — (djj — d^)/u so 
far can be tested explicitly in the first five models. For the Ising model, the Kitaev toric model, and Harper model, we have 
dt = d~ , hence a'^ = a~ ; while for both the Lipkin-Meshkov-Glick model and the Kitaev honeycomb model, and a* are 
different. 



Model 


dl 


V 


dt 


a+ 


da 


a 


ID Ising model(/ic = 1)[36] 


2 


1 


1 


1 


1 


1 


Deformed Kitaev toric model [Ac = |ln(-\/2 + l)]f68] 


In 


1 


1 


In 


1 


In 


Extended Harper model [A = 2/i]f83] 


5(or 2) 


2.5(or 1) 





2 





2 


Lipkin-Meshkov-Glick model(/ic = l)f57] 


4/3 


2/3 





2 


1 


1/2 


Kitaev honeycomb model( Jz,c = 1/2) [67, 78] 


2.50 


1 


2 


1 


2+ln 


1/2-ln 


ID AHM {tc = 0.456 for n = 2/3 )[49| 


5.3 


2.65 






1 


1.6 


Luttinger model(Ac = 1 of XXZ model) [52] 












1 



Luttinger model(Ac = -1 of XXZ model) [52] - - - 1 



tonian can be transformed into, 



H 



-XN+[2X~j~ l]a^a- 



1-7 



(at2 + . (155) 



Obviously, the Hamiltonian is quadratic, and can be di- 
agonalized via the standard Bogohubov transformation, 



= cos 

= smh{0)b'' + cosh{e)b. 



(156) 

Here b{b^) is bosonic annihilation (creation) operator, 
and the hyperbolic functions are set to ensure the bosonic 
commutation relation [b,b^ = 1. At the condition 



tanh[26'(/i > 1)] = 



1-7 



2A-7 - 1 

the Hamiltonian becomes a diagonal form 



H 



-X{N + 1) + 2v/(A-l)(A-7) (^b^ + i j 



(157) 



(158) 



Therefore, the set of eigenstate of the Hamiltonian is sim- 
ply denoted as {\n)} with 



-\{N + 1) + 2 V(A - 1)(A - 7)(n + 1/2), (159) 



where n is the number of quasiparticles. The driving 
term in the Hamiltonian becomes 

- ^ CT* = -2S, = 2a^a - 2S, 

i 

= -2S + 2[cosh(0)fe^ + sinh(6l)&] 

X [sinh(6')6t + cosh(6l)&] , (160) 

in which only sinh(29)b'^^ is the relevant term acting on 
the ground state and on projecting excited state. The 
fidelity susceptibility can be calculated, to the leading 
order, as 



XF (7,A) 



(1-7)= 



32(l-A)2(A-7)= 



(161) 



So the critical a"*" — 2. In the region < A < 1, the 
Lipkin-Meshkov-Glick model can be diagonalized simi- 
larly. The fidelity susceptibility becomes, to the leading 
order. 



Xf (7,A) 



N 



4v/(1-A2)(1-7) 



(162) 



The critical exponent a~ = 1/2. Therefore, the critical 
exponents are different on both sides of the critical point. 
The exact results are the same as those obtained from the 
numerical data. 

As a brief summary of the scaling and universality of 
the fidelity susceptibility, we collect the critical exponents 
of the fidelity susceptibility in various quantum phase 
transitions, and show them in Table. HI From the ta- 
ble, we can see that quantum phase transtions can be 
classified into two distinct types, depending on whether 
the adiabatic dimension changes or not during the tran- 
sition. For the former case, the change in the adiabatic 
dimension plays naturally a role of order parameter. 



E. Higher order of the fidelity 

The fidelity susceptibility denotes only the leading 
term of the fidelity. In case that the fidelity suscepti- 
bility is invalid, it might be useful to look into the higher 
oder term in the fidelity. In this subsection, we present 
some basic formulism of the fidely expansion. 

The overlap between two wavefunction |^'o(A)) and 
|5'o(A + 6X)) can be expanded to an arbitrary order, i.e. 



fiX,X + 5X)^l + J2 



(sxy 



4'o(A) 



Qn 



«'o(A) 
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Therefore, the fidehty becomes 



(sxr 



n=l 



oo 

E 



(5A)" / d 



n! \ dX 

((5A)'"+" 



7 ' r?, ' 



d 



(163) 



9A« / \ 9A" 



V. FIDELITY PER SITE, MIXED STATE 
FIDELITY, AND RELATED 

Mathematically, the fidehty can be roughly classfied 
into two types, the pure-state fidehty and the mixed- 
state fidelity. In this section, we introduce various fidelity 
measure under different physical conditions rather than 
mathematical view. We wiU see below, various fidehty 
expressions can be traced back to the original fidelity 
definition by Uhlmann. 



Using the relation for a given n 



E 



— ' m!(n — to)! \ 9A' 

m=o ^ ' ^ 



Qn- 



— *o = (164) 



we can find that 



= 1 - Y.(^X) 



(165) 



1=1 



where 



l—7n-\-n 



1 / d 



!n! \ aA™ 



P 



where P is defined in Eq. ([55]) . It is easy to check that 

X'p^ is zero and x^f^ the fidelity susceptibility we dis- 
cussed in previous subsections. 

On the other hand, according the perturbation theory, 
the ground-state wavefunction, up to the second order, 
is 



— £/0 — -C/n 



{Eq — Em){Eo — En) 



{Eq — E„ 



E 



2 ^0 - ^")' 



(167) 



Therefore, 



Xf 



( frOO 1 TTOn\ ttQii tttiO 



jjQii ^nm ^rnO 



{E^ - E^){E^ ~ En?- 



(168) 



Eqs. (fT66)) and p68)) conclude the main formulism of 
the higher order expansion of the fidelity. Up to now, the 
physical meaning of the high order term in the fidelity is 
still not clear. 



A. Fidelity per site 

In quantum many-body systems, the fidelity between 
two ground states increases as the system size increases. 
It usually scales like in the large L limit, where 

N — L"^ is the system size and !F is called the scaling 
parameter. Therefore, Zhou, Zhao, and Li [46l | proposed 
the scaling parameter (also called the fidelity per site) 
might be a good quantity to describe quantum phase 
transitions. Precisely, the fidelity per site is defined as 



T(X, A') = lim F 

N 



1/N 



(A, A'). 



(169) 



The expression can also be written in terms of logarith- 
mic fidelity 

lnJ^(A,A')= lim i-lnF(A,A'). (170) 

>CXD I\ 

Clearly, the fidelity per site has following properties: 

1) It is symmetric under interchange A < — *■ A', 

2) ^(A,A) = 1, 

3) < J^(A,A') < 1. 

Taking the one-dimensional transverse-field Ising 
model as an example, the fidelity between two ground 
states at h and h' is 

Fih,h') = \{^o{h')\^o{h))\ = n^o«(^'= -^5^)' (171) 

fc>0 

then the logarithmic fidelity becomes 

lnJ^ih,h') - lim l-J2cos{e,-e',), (172) 



1 



N 

fc>0 

r 

cos(6'fc - 9k)dk. 



(173) 



Therefore, one can easily discuss the critical behavior of 
the fidelity per site or the logarithmic fidelity in quantum 
phase transitions. 

The fidelity per site can be mapped onto the partition 
function of a classical statistical vertex lattice model with 
the same lattice geometry and dimension [ssj . The map- 
ping is due to the recent remarkable finding that any 
state of a quantum lattice sys tem i nay be re presented 
in terms of a tensor network Il47l . Il48j |. such as 

a matrix product state for one-dimensional systems or 
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a projected entangled-pair state for systems in D > 2 
dimensions. Then the fidelity between any two groud 
states can be cal culat ed in the context of the tensor net- 
work algorithms [1495 • Clearly, such an approach can be 
generalized to the fidelity susceptibility. 

On the other hand, in case that the adiabatic dimen- 
sion and system's dimension are not the same da ^ d in a 
certain phase, the logarithmic fidelity should be redefined 
[9^ . in principle, 

lnJ^(A,A')= lim -Lini^(A,A') (174) 

L ►oo _L " 

which keeps intensive inside the phase. As a simple ex- 
ample, in the Lipkin-Meshkov-Glick model, if A, A' < 1, 
then lni^(A, AM oc N; while if A, A' > 1, lnF(A, A') itself 
is intensive 19211. 



B. Reduced fidelity 

The reduced fidelity is defined as the fidelity between 
two reduced-density matraces corresponding to the states 
of a part of the system separated by a small distance in 
the parameter space. Precisely, if we divide the system 
into two parts: A and B, the reduced state of part A can 
be obtained by tracing out the degree of freedom of part 
B, i.e. 



PA{h)^tYB{\^oW){^o{X)\)- 



(175) 



Then the fidelity between two reduced states paW and 
PAi^) at A and A' is simply the mixed-state fidelity 

^(A, A') = try'[p^(A)]i/2p^(A')[p^(A)]i/2. (176) 

In case that the ground state is A^-fold degenerate, i.e. 
\'^}j{h)),i = 1, ..A^, one can either define the thermal 
ground state as 



pih) 



(177) 



or choose either of them as a physical ground state. The 
reduced-density matrix can be obtained in a similar way. 

The reduced fidelity approach to quantum phase tran- 
sitions was firstly discussed by Zhou 45] in order to study 
the relation between the fidelity per site and the renor- 
malization group flows. The reduced fidelity was used 
later to study the quantum phase transitions in a super- 
conductin g la ttice with a magnetic impurity inserted at 
its center [53] • Nevertheless, for a pedagogical purpose, 
here we introduce its application to another class of quan- 
tum phase transitions induced by a sequence of ground- 
state level-crossing, such as the magnetization process. 
In this case, the fidelity between two "global" ground 
states is not suitable because it drops to zero at each 
level-crossing point. For these systems, the Hamiltonian 
take typically the form 



where h is the external field and M is the magnetiza- 
tion of the ground state. Unlike the Hamiltonian of Eq. 
P3p where Hq and Hi usually do not commute with each 
other, M in Eq. (|178p commutes with Hq. Then Hq and 
M have the same eigenvectors. Suppose the magnetiza- 
tion of the system is zero in the absence of the external 
field, the system will be magnetized step by step as the 
external field increases until it is fully magnetized at a 
certain transition point he- Typically, the ground state 
has two phases, i.e., a partially magnetized phase and a 
fully magnetized phase. For the latter, the ground state 
does not change as the external field changes. Then the 
fidelity is always one and the fidelity susceptibility is zero. 
While for the former, the ground state undergoes infinite 
level-crossings and the fidelity drops to zero at each cross- 
ing point. Therefore, it is not convenient to characterize 
this type of phase transition in terms of global state fi- 
delity. The difficulty can be overcome if one consider only 
the reduced state of a local part of the system, which is 
described by a reduced-density matrix. Mathematically, 
the "orthogonality" of two pure states is very sensitive 
to the distinguishability between the s tates , such as good 
quantum numbers or which- way flag [l5Clj |. However, if 
one consider only the local part of the system, the rate of 
"orthogonality" is reduced because the reduced state is a 
mixed state that usually is free of the conserved quantity. 

The Lipkin-Meshkov-Glick model still provides a very 
good example [ll, [t^] ■ If 7 = 1 , the Hamiltonian of Eq. 
()150p commutes with the 2-component of the total spins. 
Then S"^ is a conserved quantity. The ground state can 
be written in US', S"^)}. For S = N/2, the eigenenergies 
are 

E (Af, /,) = I (^Af - - ^ (1 + h') , (179) 

and the ground state (with quantum number Mq) is then 
obtained 
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,0 < /i< 1 



(180) 



where I{x) gives the integer part of x. We can define the 
reduced state of a single spin as 



P = 



1 / 1 











1 



(181) 



The fidelity between two reduced states around the each 
crossing point is 



H{h) = Hq~ hM, 



(178) 



+ 2v/(l-(^.>^)(l-(^.>^'+^), (182) 

where (uz)"' denotes the expectation value of az of jth 
state during the level-crsossing process. The fidelity sus- 
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ceptibility can be defined in the similar way, in the ther- 
modynamic hmit, 



-21nF 



(183) 



where 5h = hj+i — hj is the difi^erence of h between two 
level-crossing points. 



Similarly, if the driving term in the Hamiltonian is a 
Zeemann-like term, say —hM, which is crucial in Landau- 
Ginzburg- Wilson symmetry-breaking theory, then the fi- 
delity susceptibility is simply the magnetic susceptibility 

xSi, 
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(193) 
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C. Fidelity between thermal states 

At finite temperatures, an equilibrium state of a ther- 
mal system is described by a mixed state [l5ll |. 



z ^ 



(184) 



where /? — 1/T is the inverse temperature with Boltz- 
mann constant kg = 1, Z is the partition fmiction de- 
fined as 



(185) 



The thermal fidelity susceptibility is similar to the 
ground-state fidelity susceptibility. Both of them are a 
kind of structure factor because of 
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(194) 
(195) 



in thermodynamics. The main difference is that the 
ground-state fidelity susceptibility involves the dynamic 
behavior of the system. 

As a simple application, we take the two-dimensional 
Ising model on a square lattice as an example. The 
Hamiltonian reads, 



and l^^n) is the eigenstate of the system's Hamiltonian 
with eigenvalue En- Therefore, if we choose the temper- 
ature as a driving parameter, the fidelity between two 
thermal states at /3 - (5/3/2 and /3 5/3/2 is [13] 



F{P- 6(3/2, (3 + 6/3/2) 
tr^p{P-60/2)p{p + 6/3/2) 
Z{[3) 

^Z{(3~5p/2)Z{/3 + 5(3/2)' 



(186) 
(187) 



This result is very interesting. It establishes a relation 
between the state overlap and well known thermal quan- 
tities. Therefore, one can understand the state evolution 
at finite temperature from the knowledge of thermody- 
namics. Since the Helmholtz free energy is 
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the thermal fidelity susceptibility then becomes [43, |4; 
-21nF 



XF 



[2Gi(3) - G(/3 + 6(3/2) - G{(3 - 5/3/2)] 
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To obtain the above result, we have used the following 
standard relations in the thermodynamics 151 



dG_ 
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where the sum is over all pairs of nearest-neighbor sites 
i and j, and the coupling is set to unit for simplicity. 
The model is certainly the most thoroughly researched 
model in statistical physics [206l i207i | . The results for a 
40 X 40-site system are shown in Fig. [121 Clearly, there is 
a maximum point in the line of the specific heat, whose 
scaling behavior to an infinite system defines the critical 
point. Meanwhile, the middle picture in Fig. [T^j shows 
various fidelity calculated from different temperature in- 
terval. This obvious difference in the fidelity disappears 
if we distill the fidelity susceptibility from them, as shown 
in the right picture of Fig. [T^ 



D. Operator fidelity 

The operator fidelity was proposed by Wang, Sun and 
Wang [6J|. For two arbitrary linear operators A, B de- 
fined in a d-dimensional Hilbert space, the expectation 
value of their product are tr(^i?), which is a kind of in- 
ner product in the d^-dimensional (or cP — \ for Hermitian 
operators) space. On the other hand, in a d dimensional 
Hilbert space, any state acted after a linear operator be- 
comes another state in the Hilbert space. Then the fi- 
delity between states can be generalized to the operator 
level. Specifically, for two unitary evolution operators Uq 
and C/i, the fidelity between them can be calculated as 
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tr{UlU,)\. 



(197) 
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FIG. 12: (Color online) LEFT: The specific heat as a function of the temperature T for 40 x 40 Ising model. MIDDLE: The 
fidelity between two states separated by different temperature interval 5T = 0.02,0.04,0.06.0.08,0.1,0.15,0.20 for lines from 
up to bottom. RIGHT: the fidelity susceptibility xf as a function of T, obtained from the data of the middle picture. All lines 
in the middle picture collapse onto a single line (From Ref. [i^). 



Precisely, for a general Hamiltonian, the evolution oper- 
ator can be expressed as 

U(t) = l-iSX [ dtiV{ti) 
Jo 

-{5Xf [ dti r dt2V{ti)V{t2) + --\m) 
Jo Jo 

where V{t) = exp{iHQt)V exp{—iHot) denotes the per- 
turbation term in the interaction picture. The trace of 
the evolution operator is, 

t^[U{t)] = l-iSXElWit)]- (199) 
-^-^-r[W{tf]+---, (200) 

where 

tF(A) = itr(A), (201) 

and 

W{t) = [ dtiV{ti). (202) 
Jo 

The left hand side of Eq. p99p denotes the inner product 
of two states at different time. Then the fidelity becomes 

F^ = l- {SXf [tF [W{tf] - (MWit)f] +■■■ . (203) 

Similarly, the operator fidelity susceptibility can be ex- 
tracted 

XF = [W{t)^] - [t^W{t)] ^ . (204) 



The operator fidelity is remarkable. In its approach to 
quantum phase transitions, it is a good indicator regard- 
less of the ground-state degeneracy. Moreover, it also re- 
veals that in the state evolution, the driving mechanism 
is due to the fluctuation of W{t). 



E. Density-functional fidelity 

The density functional theory (DFT ) [208l . [209} is 
based on the Hohenberg-Kohn theorem j208| , which as- 
serts that the ground-state properties are uniquely de- 
termined by the density distribution that minimizes 
the functional for the ground-state energy EqIux]. To 
date, the DFT becomes the most successful method for 
first-principle calculations of the electronic properties of 
solids. Since the normalized density distribution captures 
the most relevant information about the ground state, Gu 
[s^ l tried to link the fidelity and quantum phase transi- 
tions via the DFT. 

For a general Hamiltonian system 

HiX)=Ha + \HI + Y,^^xn,, (205) 

X 

where Hj is the interaction term and A denotes strength, 
and fix is the local (pseudo)potential associated with 
density distribution {ux}. The index x can be discrete 
or continuous depending on the system under study. 
Though in the local-density-approximation (LDA) calcu- 
lation, Tlx usually refers to the density of electrons in real 
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space, it can also be generalized to the population in con- 
figuration space of a reduced-density matrix or the den- 
sity of state in energy(monientum) space. The density 
distribution can be obtained by the Hellmann-Feymann 
theorem 



*o(A) 



dH 



*o(A) 



(206) 



where 1^*0 (A)) is the ground state, and J2x = ^- The 
density-functional fidelity is defined as the distance be- 
tween two density distributions at A and A' in the pa- 
rameter space 
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F{\ A/) ^ trVn(A)n(A'). 



(207) 



Since the density distribution is experimentally measure- 
able, the density-functional fidelity then provides a prac- 
ticable approach for experimentalist to study quantum 
phase transitions in perspective of information theory. 

Expanding the density-functional fidelity to the lead- 
ing order, one can find the density- functional fideltiy sus- 
ceptibility has the form 



XF 



E- 



(208) 



FIG. 13: (Color online) The fidelity of the one-dimensional 
transverse-field XY model as a function of h and 7, for a 
system of A'' = 100 and SX = S'y — 0.1. The colored curves on 
the F — plane constitutes a contour map. The deep grooves 
in the curved surface of the fidelity separate naturally the 
whole region into three phases: one order phase in the middle 
and two polarized phases on both sides (Reproduced from 
Ref. Hi]) 



A. Pure-state fidelity 



One- dimensional spin systems 



Therefore, if one regards dn^/dX as an independent func- 
tion besides the density distribution n^, the density- 
functional fidelity susceptibility is a functional of Ux 
and dux/dX, both of which, in principle, maximize the 
density-functional fidelity susceptibility at the critical 
point. 

Mathematically, the density-functional fidelity and 
thermal-state fidelity is related to the Fisher-Rao dis- 
tance. For the probability {hx} defined in parameter 
space {Aa}, the Fisher- Rao distance is defined as 



dnxdn-. 



The Fisher-Rao metric then can be written as 
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1 dux dux 
Ux dXa dXb 
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VI. 



FIDELITY IN STRONGLY CORRELATION 
SYSTEMS 



The fidelity approach to quantum phase transitions has 
been applied to many strongly correlated systems besides 
the one-dimensional transverse-field Ising model and the 
Lipkin-Meshkov-Glick model. This section is devoted to 
a survey of fidelity in these systems. We mainly focus 
on those models which are well studied and whose phase 
diagrams are known. 



The one- dimensional transverse-field XY model: The 
model is an extended version of the one-dimensional 
transverse-field I sing model. The XY model can be ex- 
actly solved too [l29l Il30l | . Its Hamiltonian reads 



^(7, A) = -J2 



1 + 7 



1-7 



V V 



Xa: 



(211) 



where 7 defines the anisotropy and A represents the ex- 
ternal magnetic field. Obviously, if 7 = 1, the XY model 
becomes the transverse-field Ising model which is pre- 
sented as an example in section Hill The role of quantum 
entanglement in the quantum phase transitions occurred 
in the XY model has been widely studied 0, H, fl^Jlj | . 
It was shown that either the pairw ise entang lement [7|,|8|, 
as measured by the concurrence |l52l Il53l | or the two- 
site local entanglement [13, [IH shows interesting singu- 
lar and scaling behavior around the critical point. The 
fidelity approach to the one-dimensional transverse-field 
XY model was firstly done by Zanardi and Paunkovic 
[3^ . They used the XY model as an example to present 
a new characterization of quantum phase transitions in 
terms of the fidelity. Their physical intuition behind the 
fidelity approach is in order: since the fidelity measures 
the similarity between two states, then the fidelity be- 
tween two ground states obtained for two different values 
of external parameters should has a minimum around the 
critical point (See Fig. [T5)l . The work is one of original 
works in this promising field. 
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FIG. 14; (Color online) The theoretical (line) and experimen- 
tal (asterisks) fidelity in the quantum phase transition of the 
Ising dimer induced by the ground-state level-crossing. (From 
Ref. [z3]) 
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FIG. 15: (Color online) The experimental fidelity of the Ising 
dimer measured by the nuclear-magnetic-resonance experi- 
ments for Sh^ = 0.2 (*) and Sh^ = 0.3 (x). (From Ref. 
0) 



Meanwhile, Zhou, Zhao, and Li [4^ studied the critical 
phenomena in the XY model in terms of the fidelity per 
site. They found that the logarithmic function of fidelity 
per site with respect to the transverse field is logarith- 
mically divergent at the critical point. The scaling and 
universality hypothesis based on the fidelity per site was 
also confirmed in the vicinity of the transition point. This 
work is the first one that proposed the idea of fidelity per 
site. 

A noticeable achievement is the experimental detec- 
tion of the quantum phase transitions in terms of the 
fidelity [7l|. Using the nuclear-magnetic-resonance tech- 
nique, Zhang et al measured the sensibility of the ground 
state of the Ising model to perturbations when it comes 
to the critical point. Their system can be described by 
the Hamiltonian 

= alal + B.,{al + ) + B,(al + al). (212) 

If Bx = 0, the system undergoes a first-order phase tran- 
sition induced by the ground-state level-crossing. They 
found that the ground-state overlap shows a drop at 
the crossing point (|14p . While if = 0, the model 
becomes the transverse-field Ising model. Though the 
second-order quantum phase transition occurs only in the 
thermodynamic limit, the fidelity in a small sample can 



still tell us the significant change in the structure of the 
ground state around the critical point (See Fig. [T5)) . 

The one- dimensional XXZ model: The Hamiltonian of 
the one-dimensional XXZ model reads 

L 

H{X) = (^^^^+1 + "'"j+i + K^l+i) • (213) 

The XX Z model can be solved by the Bethe-ansatz 
method |l54l Il55l Il56l | , through which the energy spec- 
tra can be fully determined. The ground-state of the one- 
dimensional XXZ model consists three phases. If A < — 1, 
all spins are align to the same direction, the ground-state 
is in a fully polarized state; when — 1 < A < 1, the quan- 
tum fluctuation term dominates, then the ground state is 
a quantum fluctuation phase; while if A > 1, the antifer- 
romagnetic coupling dominates, the ground state is in an 
antiferromagnetic state. So there are two critical points 
A = ±1. There are various tools to witness the quantum 
phase transition occurred at A = 1. They include the van- 
ishing energy gap and divergent correlation length Il57[ 
as A ^ l"*", abrupt change in the spin stiffness [l58lll5g |. 
and the maximum value of quantum entanglement |l2j|. 

The fidelity approach to the quantum phase transition 
of the one-dimensional XXZ model is not easy because 
the ground-state wave function is not known. Yang (s^ 
first used the Luttinger Liquid model to describe the one- 
dimensional XXZ model in the quantum fluctuation re- 
gion [l60j |. i.e. 



iJ(A) = ^ J dx (^KU{xr + ^idx<i>f 



Here 
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and H, $ are bosonic phase field operators. They ob- 
tained the fidelity as 



F(K,K') = TT , , 

and the fidelity susceptibility 

Xf(A) 1 1 



L 



4:[k - arccos(A)]2 1 - A^ ' 



(217) 



(218) 



Therefore the critical exponent a ~ I. So they conclude 
that the fidelity might be able to signal the Beresinskii- 
Kosterlitz-Thouless phase transition occurring in the 
XXZ model. This conclusion is a little surprising. Be- 
cause in previous studies, some groups claimed that the 
Beresinskii-Kostcrlitz-Thoulcss transitions cannot be sig- 
nalled from the singularity of the fldelity susceptibility 

MM- 
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FIG. 16: (Color online) The fidelity and fidelity susceptibility 
in the spin-one anisotropic model (From Ref. ffd^]. 



The spin-one anisotropic modetThe Hamiltonian of 
the spin-one anisotropic model reads 



HiX) 



E 



(219) 

where S^^k = x,y,z) stands for spin-1 operators, and 
A and D parameterize the Ising-hke and the uniaxial in- 
teraction. The ground -state pha se diagram of the model 
consists of six phases |l 6 ll Il62l |. Specifically, for three 
cases of A = 2.59,1, and 0.5, the system undergoes 
second-order, third-order, and fifth-order quantum phase 
transitions, respectively. For the case of A = 2.59, the au- 
thors found the fidelity susceptibility as a function of D 
around D ~ 2.30. For the second case of A = 1, though 
the second-order derivative of the ground-state energy 
does not show singular behavior around D = 0.95, the 
fidelity susceptibility is still able to signal the transition. 
For the third case of D = 0.5, however, both the fidelity 
susceptibility and the second derivative of the ground- 
state energy is continuous in the critical region. The au- 
thors then conclude that the fidelity susceptibility might 
not be able to characterize the Gaussian transition. (Sec 
Fig. [TT 



The spin-ladder model: The Hamiltonian of a general 



H — y~^[J(.S'i jS*! j + i + »5'2j.S'2j-|-l) + JrSljS2, 

+y{Si,jSi^j+i){S2,jS2,j+i) 



+ Jd{SljS2J + l + S2,jSij + i) 

+K[{Si^jS2,j-\-i){S2jSi^j+i) 

-iSl,jS2,j)iSij + iS2,j + l)], 



(220) 



where the indices 1 and 2 distinguish the lower and upper 
legs of the ladder and i labels the rungs. The ground 
state of the spin ladder is very complicated. However, if 
we redefine the parameters 
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—u, K = Jj. 
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Jd = 0, (221) 
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the Hamiltonian becomes a one-parameter model. The 
ground state of the system can be explicitly written as a 
matrix product state 

l^o(u,w')> = -^=ir[gi{u)g2{u') ■ ■ ■ g2N-i{u)g2N{u')], 



/N, 



where Nc is the normalization constant. Such a ground 
state undergoes two second-order quantum phase tran- 
sitions at u = and u = oo. At u = 0, the ground 
state changes from the dimerized phase to the Haldane 
phase. The latter can be described by an effective Hamil- 
tonian of the S = 1 Affleck-Kennedy-Lieb-Tasaki chain 
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FIG. 18: (Color online) The dependence of the fidelity sus- 
ceptibility of the Hubbard model on the system size at f/ = 
(From Ref. Es 



|164| . Fig. [T7] shows the fidelity approach to the phase 
transition occurring at u = 0. Clearly the fidelity shows 
a minimum and the fidelity susceptibility shows a sharp 
peak at M = 0, corresponding to an abrupt change in the 
ground-state wavefunction. 



2. One-dimensional fermionic systems 

The one- dimensional Hubbard model: Th e H amilto- 
nian of the one-dimensional Hubbard model [l65l | reads 

L L 

HnM = -Y^ ^t4,c,+5,, + C/^n,-|nj-x, (224) 



j=l 5=±1 



where cj ^ and Cj,a, cr =T7 i ^-re creation and annihilation 
operators for fermionic atoms with spin a at site j respec- 



tively, 71c 



J, and U denotes the strength of on-site 



interaction. Besides the obvious SU(2) symmetry i n th e 
spin sector, the model has c harge SU(2) symmetry |l66l |. 
The global symmetry jl67l | of the model is SO (4) since 
half of the irreducible representations of SU(2)(8)SU(2) 
are excluded. The H ubbard inodel was solved by the 
Bethe-ansatz method [l6l HH, UzO]. For the half-filled 
case, the system undergoes a Mott-insulator transition 
at the critical point Uc = 0- According to the exact solu- 
tion, the ground-state energy can be differentiated to an 
arbitrary order, therefore the transition is of infinite or- 
der [l69l. [T70|(Beresinskii-Kosterhtz-Thouless like). The 
role of quantum entanglement in the phase transitions oc- 
curred in the Hubbard model was addressed by Gu et al 
[13], Deng et al, [lB| and Larsson and Johannessonp^. 
The fidelity approach to the one-dimensional Hubbard 
was firstly studied by You et al [i^. They found that 
the fidelity susceptibility is not singular at the transition 
point and concluded that it might not be able to signal 
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FIG. 19: (Color online) A sketch of the asymmetric Hubbard 
model which can be realized in an optical lattice. The solid 
sinusoid denotes the periodic potential formed by two interfer- 
ing standing laser waves. Two species of fermionic atoms are 




FIG. 20: (Color onhne) The scaling behavior of the FS as 
a function of t for the cases of n = 2/3 (LEFT) and n = 1 
(RIGHT). Here U = 30. (From Ref. [4|) 



the transition. Especially at the critical point Uc = 0, 
the Hamiltonian is diagonal, and the fidelity susceptibil- 
ity per system length, as an intensive quantity, is a con- 
stant (See Fig. [T5|) . The fidelity in the Hubbard model 
was later revisited by Venuti et al f63|, with the strat- 
egy of bosonization and Bethe-ansatz techniques. They 
showed that the metal-insulator phase transition can be 
insightfully analyzed in term of the fidelity. The fidelity 
susceptibility shows divergences depending on the path 
approaching to the critical point. Bosonization results 
shows that the fidelity susceptibility may diverge as 

XF oc C/"", 

if the doping rate is proportional to -v/C/exp]— 27r/[/]. 
The authors also performed exact diagonalization for a 
system up to 14 sites. Based on the scaling analysis, their 
results indicate that the fidelity might be super-extensive, 
hence divergent in the thermodynamic limit. However, 
they also stated that the available sizes of the exact diag- 
onalization are too small to provide a conclusive answer. 
Therefore, whether the fidelity susceptibility shows sin- 
gular behavior in the one-dimensional Hubbard model is 
still not conclusively answered. 

The one- dimensional asymmetric Hubbard m odel: The 
Hamiltonian of the asymmetric Hubbard model |l7ll . ll72| 
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reads 



j=l S=±l a 



(225) 

where is cr-dependent hoping integral. Not much at- 
tention was paid to the model in the last century be- 
cause we did not have a realistic system that the model 
can be applied. However, in describing a mixture of two 
species of fermionic atoms in optical lattices which has 
been realized by recent experiments on the cold atoms 



173 1 , the model itself becomes a current research interest 
174 Il75l Il76j . The ground-state phase diagram of the 



asymmetric Hubbard model can be understood from its 
two lin iiting cas es, i.e. the Falicov-Kimball (FK) model 
region 1 1771 . Il78| {ti = 0) and the Hubbard model region 
(HM) [l65j | (if = ti). The schematic phase diagram of 
the mod el wa s obtained from the renormalization group 
analysis [174| . Subsequently, a quantitative phase di- 
agram was also obtained by analyzing structure factor 
with the density- matrix renor malization group 175] and 
bosonization techniques [l76j |. The fidelit y a pproach to 
the model was firstly done by Gu et al [4^. The au- 
thors found that the fidelity susceptibility can be used 
to identify the universality class of the quantum phase 
transitions in this model. That is the quantum phase 
transitions occurred at different band filling can be char- 
acterized by the critical exponents of the fidelity suscep- 
tibility. Fig. [20] shows the fidelity susceptibility as a 
function of ij,/^]- for the cases of n = 2/3 and n = 1, 
and various system sizes. Obviously, the fidelity suscep- 
tibility diverges quickly as the system is away from the 
half-filling, and relatively slow at half-filling. The scaling 
analysis reveals that the maximum value of the fidelity 
susceptibility scales like 
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at n = 2/3 case, while 
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at half-filling. Then the intensive fidelity susceptibility 
shows singularity away from half-filling, while no singu- 
larity at half-filling. These observations support their 
conclusions on the role of fidelity in describing the uni- 
versality class. 

The one- dimensional Bose-Fermi mixture: The sim- 
plest one-dimensional Bose-Fermi Hubbard model can be 
modeled by 

N 

Hbf = -J^i^F^c^+i+tBblbi+i+R.c.) 

i=l 

N N 

+ubf J2 + ubb J2 - m^) 



Here hi {h\) and Ci (c|) are the bosonic and fermionic 
annihilation (creation) operators at site i, respectively. 




FIG. 21: (Color online) The ground-state phase diagram de- 
fined on the Ubb-Ubf plane of the one-dimensional Bose- 
Fermi Hubbard model in terms of the fidelity (Reproduced 
from Ref. [H). 



^F(iB) is the hoping integral of fermions (bosons). Ubf 
denotes the on-site interaction between fermion and 
bosons, and Ubb the interaction between bosons. The 
ground state of the model has been studied by the quan- 
tum Monte-Carlo method J 798] . Several phases, includ- 
ing Luttinger liquid phase, density wave phase, phase 
separation state, and Ising phase, are predicted. The 
fidelity approach to the quantum phase transitions oc- 
curring in the ground state was done by Ning et al [s^ ]. 
As shown in Fig. [^TJ Luttinger liquid phase, density wave 
phase, phase separation state can be observed from the 
behaviors of the fidelity. However, it is difficult to find 
the phase boundary between the Ising phase and density 
wave phase. The authors interpret that the phase tran- 
sition between the Ising phase and density wave phase 
is of the Beresinskii-Kosterlitz-Thouless type since in the 
large Ubb and Ubf limit, the effective Hamiltonian of 
Eq. (|226p becomes the one-dimensional XXZ model. In- 
stead, they found that the concurrence 152, 153], as a 
measure of entanglement between two 1/2 spins, can lo- 
cate the transition point because the concurrence shows 
a maximum at the isotropic point of the XXZ model [l3] . 
Nevertheless, as we mentioned before, the role of fidelity 
in the Beresinskii-Kosterlitz-Thouless phase transitions 
is still controversial. 

The ex tended Harper model: The Harper model jlSOl . 
llSlL Il82| ] was proposed to describe electons in a two- 
dimensional periodic potential under a magnetic field, 
for u nderstanding Hofstadter-butterfly energy spectrum 
[l83j ] . For a system of electrons moving on a triangle lat- 
ticein a magnetic field, the Hamiltonian can be written 
as[l82j 



Harp 



n 

+ i^g2Ti0(n+l/2)-jfc^ 

n 



(227) 
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is shown in Fig. [22lja) jl82{ . Therefore, in addition to 
metal-insulator transitions, there is also an interesting 
metal-metal transition. 
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The fidelity approach to quantum phase transitions be- 
tween the three phases of the Harper model was done by 
Gong and Tong f83|. They studied the fidelity between 
the ground state at (A,/i) and various reference states. 
For example. Fig. [22y b) shows the fidelity of a reference 
state at (Aq, /io) — (2.0, 1.0) which is a tricritical point in 
the phase diagram. At that point, the fidelity is 1, while 
away from the tricritical point, the ground state changes 
quickly, then the fidelity show distinct behavior in the 
parameter space. One can find that the ground-state 
phase diagram can be sketched out by the contour map 
of the fidelity. Such a property also is consistent with 
the primary motivation of the fidelity approach. The 
authors studied also the fidelity susceptibility along var- 
ious evolution directions. Fig. ISST c) shows the fidelity 
susceptibility in the \ — \i plane along (l/-\/5, — 2/-\/5) 
direction. In this case, the fidelity susceptibility are a 
combination of two elements in the quantum geometric 
tensor. Since the direction is not parallel to any transi- 
tion line, the fidelity susceptibility reaches maximum at 
the boundary between three phases. Interestingly, the 
authors found the critical exponents of the fidelity sus- 
ceptibility depend on the different choice of system size. 
For example, if the system size equal to Fibonacci num- 
ber Fm{F„i = Fm-1 + Fm-2) with m — 31 + 1, the adi- 
abatic dimension at the critical point = 4.9371 and 
ly = 2.4718; while form m ^ 31 + 1, ^ 2 and v = 1. 
Therefore, the critical exponent a — d^^al^ = 2 in both 
cases. 



3. The fidelity in topological quantum phase transitions 



FIG. 22: (Color online) (a) The ground-state phase diagram 
in the A— plane. (b)The fidelity F{X, jj,; Ao, /lo) for (Ao, /io) ~ 
(2.0, 1.0). (c) The fidelity susceptibility and its contour map 
as a function of A and fi along the direction (l/\/5, — 2/\/5). 
(From Ref. HI 



where ta,th,tc are hoping amplitude for each bond on 
the triangular lattice, n lattice index, 0/2 magnetic fiux 
piercing each triangle, and ky the momentum in y di- 
rection. Since there is no interaction between elec- 
trons, the model can be solved exactly. The ground 
state of the model consists of three different phase, two 
conducting phases and one insulating phase. Defining 
A = 2tb/ia,/i = tc/ta, the ground-state phase diagram 



Topological phase transitions are very special in quan- 
tum critical phenomena. They do not rely on any local 
order parameters nor on a symmetry breaking mecha- 
nism, hence cannot be described by Landau-Ginzburg- 
Wilson paradigm. A typical example of these novel 
phases is the fractional quantum Hall state, in which elec- 
trons in two dimension are strongly correlated and their 
fiuctuations are entirely quantum in nature, therefore, 
the Landau- Ginzburg- Wilson theory, which is based on 
a classical local order, might fail. The fidelity approach 
provides an alternative method to study these fascinating 
phase transitions. 

The deformed Kita ev toric model: The deformed Ki- 
taev toric model [184| is defined on a square lattice with 
spin- 1/2 degrees of freedom residing on the bonds (Fig. 



30 



t» — 




FIG. 23: (Color online) A sketch of the Kitaev toric model: 
The square lattice (LEFT) is defined on a torus. All 1/2 spins 
locate at the middle point (solid dot) of bonds. Four neigh- 
boring spins interact with other by plaquette or star operators 
(RIGHT) depending on their locations. 



FIG. 24: (Color online) A sketch of the Kitaev honeycomb 
model: Spins locate at the vertices of a honeycomb lattice 
(LEFT). Each spin interacts with three neighboring spins 
through three types of bonds, i.e. "x{y, z)" bonds depend- 
ing on their direction (RIGHT). 



left). The Hamiltonian reads 
Hdkt = — Ao ^ Bp — Ai ^ As 



Hamma, and Zanardi [66 
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Here Ag — Hies '^f ^^'^ — Yliep '^f ^^'^ ^^'^ 
plaquette operators (Fig. [221 right) of the Kitaev model, 
Ao,i > and is a parameter tuning the system across 
a topological phase transition. The ground state of the 
defor med Kitaev toric model has been obtained exactly 

\Mm - (230) 



Z{p) = ^exp 
geG 



,eG vm 
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where G is the Abelian group of all spin-flip operators 
obtained as products of star-type operators. The fidelity 
approach to the topological phase transition occurring in 
this model has been done [gl, HI] . 

The ground state of the model consists of two distinct 
phases [185i] . If /3 = 0, the model is pure Kitaev toric 
model and its ground state is a closed string condensed 
phase, in which each x{z) string is a collection of spins 
that are flipped in the (j^{a^) basis. While if (3 is very 
large, the ground state favors a fully polarized phase. 
A quantum phase transition is found to be occurred at 
the critical point /3c = (l/2)ln(\/2 + 1). The system 
shows a topological order \i (3 < (3c- The fldelity between 
two states at the points P ztSf3 was obtained by Abasto, 



F{P- 6(3/2, (3 + 513/2) 
{^oW-S(3/2)\M/3 + S(3/2)), 
exp[-/?X:»gf(g)] 



5(3 /2)Z{(3 + 5(3/2) 



(232) 
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and the fidelity susceptibility is 



-21ni^ 



XF 



Cy 



5/3^0 



(234) 
(235) 



where Cv denotes the specific heat of the two-dimensional 
Ising model. Then the fidelity susceptibility shows a loga- 
rithmic divergence at the critical point (3c. These findings 
are very interesting. The corresponding to the thermal 
phase transition occurring in the two-dimensional Ising 
model reveals that the topological phase transition in the 
deformed Kitaev toric model could be detected by the lo- 
cal magnetization 1^18^]. 

The Kitaev honeycomb model: The Kitaev honeycomb 
model was also introduced by Kitaev in search of topo- 
logical order and anyonic statistics. The model is as- 
sociated with a system of 1/2 spins which are located 
at the vertices of a honeycomb lattice (Fig. [211 left). 
Each spin interacts with three nearest neighbored spins 
through three types of bonds, called "^(y, z)-bonds" de- 
pending on t heir direction (Fig. [24l right). The model 
Hamiltonian Il84ll is as follows: 



x-bonds y-bonds 
2-bonds 

= —JxHx — JyHy — JzHz, (237) 
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FIG. 25: (Color online) The fidelity susceptibility as a func- 
tion of Jz along the dashed line shown in the triangle for 
various system sizes L = 101, 303, 909. Both up insets corre- 
spond to enlarged picture of two small portions (From Ref. 



where j, k denote two ends of the corresponding bond, 
and Ja{a = x,y,z) are coupling constants. The ground 
state of the Kitaev honeycomb model consists of two 
phases, i.e., a gapped A phase with Abelian anyon ex- 
citations and a gapless B phase with non-Abehan anyon 
excitations. A quantum phase transition occurred be- 
tween these two phase is believed to be a topological 
phase transition because no local operators can be used 
to describe such a transition. 

The Kitaev honeycomb model can be diagonalized ex- 
actly in the vortex-free subspace. The ground state can 
be written as 



V2 



Aq -I- ieq 



-q4 



|0), (238) 



where 



qx(y) 



Jx COS q,j: + Jy COS Qy + Jz, 

Jx sinqx + Jy sinqy. 

2mr L-l L-l 



(239) 
(240) 

(241) 



and a. 
sing le 



with 



-q.i(2) &re Majorana operators for two sites of a 
bound. Then the fidelity between two states is 



F^^JJcos^ (e^-e'^) 



(242) 



cos (26'q) 

COS (2^;,) 



|3.,sin(20q) = ^, 



-3-,sin(20i) 
q 



E' 



(243) 



The fidelity depends on the positions of two states in the 
parameter space. Therefore, in order to extract the fi- 
delity susceptibility, we must know the direction of line 
connecting the two points. If we define the ground-state 
phase diagram on the plane Jx + Jy + Jz = 1, and con- 
sider a certain line Jx — Jy along which the ground state 
of the system evolves at zero temperature. The fidelity 
susceptibility becomes [i^l 



XF 



16 ^ 



sm qx 



(244) 



Fig. shows the fidelity susceptibility's dependence on 
the driving parameter (red line in the figure) for vari- 
ous system size. The authors also performed the scaling 
analysis, and found that the fidelity susceptibility scaling 
hke 



XF 

L2 



1 



(J, -1/2)1/2' 



(245) 



around the critical point Jz — 0.5+ . While in the gapless 
phase, Gu and Lin (t^ found that the fidelity suscepti- 
bility is superextensive, i.e. 



XF 



L^lnL. 



(246) 



Then the critical exponent is quite different at the other 
side of the critical point, i.e. 



XF 



J.-l/2|i/2 



InL 



ln|J,-l/2|, 



(247) 



Meanwhile, the topological phase transition was stud- 
ied in terms of the fidelity per site by Zhao and Zhou [6^ . 
According to the definition of Eq. (|170p , the fidelity per 
site in the ground state of the Kitaev honeycomb model 
has the form 



\nd{J,J') = -l^ln[cos(0q-0;)] , 
q 

(27r)^ Jo Jo 



(248) 



(249) 



where 9q is determined by Eq. (|243p . The above expres- 
sion makes it possible to study the scaling and critical 
behavior of the fidelity per site easily. More precisely, if 
J' is fixed, they found that the fidelity per site is loga- 
rithmically divergent when J is varied such that a critical 
point is crossed. For example, if only Jx of J is consid- 
ered as a driving parameter, they found the fidelity per 
site 



d2 lnd(J, J') 



dJ2 



InL, 



(250) 
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(1^;J) p,:(gJ+5J) 
Cbulk;J) Pj:(bulk;J+5J) 
(1+1;J) p :(l+l;J+SJ) 



FIG. 26: (Color online) The single-site reduced fidelity in the 
two-dimensional s-wave conventional superconductor. (From 
Ref. EE 



where Jxm denotes the position of extremum point for a 
finite system, and 



\nd{J,J') 



In I Jj; - Ja;c|, 



(251) 



in the thermodynamic limit. 

Therefore, the fidelity now is beUeved to be able to 
characterize the topological phase transition occurring in 
the Kitaev honeycomb model [H, [g^I • 



B. Mixed-state fidelity 

1. Reduced fidelity in quantum phase transitions 

There are several works on the role of reduced fidelity 
in quantum phase transitions. 

The one- dimensional transverse-field XY model: The 
reduced fidelity approach to the XY model was firstly 
touched by Zhou [45| in order to confirm the relation be- 
tween the fidelity per site and renormalization group flow. 
He noted that the fidelity per site can capture nontrivial 
information including stable and unstable fixed points in 
the renormalization process in which the entanglement 
cannot. 

Ma et al [8l| and Li [s^ studied the two-site reduced 
fidelity in the ground state of the transverse-field Ising 
model independently. In the work by Ma et al, the fidelity 
of the reduced state of two neighboring site in the model 



F - tr^p(A)i/2p(A + a)p(A)i/2, 
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where p{X) is the reduced-density matrix of the two spins. 
Then the reduced fidelity susceptibility is defined by 



XF 



lim 



-21nF 



(253) 



For this model, p(A) can be expressed in terms of the 
two-site correlation functions that can be calculated ex- 
plicitly. This property enables them to study the scaling 
and critical behavior of the reduced fidelity easily. They 
found that, for a finite A^-site system, the function 

XriXm, Nf'^ - xf(A™, Nf/^ = Q[N''{\ - A„)] (254) 

is unversal around the critical point. The scaling rela- 
tion reveals that reduced fidelity susceptibility diverges 
logarithmically 



Xf{X) 



1/2 



InlA- A, 



(255) 



which is quite different from the global state fidelity in 
the same model. While Li 88] calculated the reduced 
fidelity directly. It is reported that the extremum of the 
reduced fidelity scales like F oc (InA^)^-^'^ around the crit- 
ical point. Later, You et al [90| also make an extension 
to the one-dimensional XY model. Similar results are 
obtained. 

The Lipkin-Meshkov-Glick model: The quantum phase 
transition occurring in the Lipkin-Meshkov-Glick model 
was also studied in terms of the reduced fidelity 0, [tI] . 

It is not convenient to use the global-state fidelity to 
characterize those quantum phase transitions induced the 
continuous level-crossing, such as the magnetization pro- 
cess, because the global-state fidelity drops to zero at 
each level crossing point. Kwok et al proposed that 
the reduced fidelity can overcome the difficulty. They 
used both the Lipkin-Meshkov-Glick model and the one- 
dimensional XXX model as examples to show that the 
reduced fidelity and its leading term, i.e. the reduced fi- 
delity susceptibility help to study scaling and critical be- 
havior. Ma et al 73^ also study the behaviors of the two- 
site reduced fidelity in the Lipkin-Meshkov-Glick model 
model. 

The one- dimensional extended Hubbard model: 
Though the role of the ground-state fidelity in the 
Hubbard model is still not very clear, it is shown by 
Li [il] that the two-site reduced fidelity is able to 
sketch out the ground-state phase diagram of the model. 
Previously, it is firstly reported by Gu et that the 

contour map of single-site entanglement can describe the 
change of symmetry in the ground-state of the system. 
Due to the competion betwee various correaltion, the 
entanglement usually reaches a maximum at the critical 
point; while from fidelity approach, it is shown that 
fidelity shows a minimum at the critical point due to the 
abrupt change in the structure of the ground state. 

Two-dimensional s-wave conventional superconductor: 
The Hamiltonian reads 



H 



E ^4 



Sf 



(256) 



- ^ J (cos ipcl„al„,CQ„, + sin ipcl^al^, 
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FIG. 27: (Color online) The scaling behavior the reduced 
fidelity susceptibility of the one-dimensional transverse-field 
Ising model. (From Ref. [sj] ). 

The first three terms in Eq. (|256p denote the Hamilto- 
nian of s-wave conventional superconductor, and the last 
term represents the interaction between electrons and an 
classical spin placed at the origin. Paunkovic et al [s^ 
first studied the reduced fidelity (or the partial-state fi- 
delity) in the ground state of the system. They observed 
that the one-site reduced fidelity shows a sudden drop 
in the vicinity of the quantum phase transition. Since 
the one-site reduced-density matrix depends on a single 
quantity, i.e. magnetization. They interpret the drop of 
the reduced fidelity due to that the on-site magnetization 
plays an order parameter in this model. Such a behavior 
is very similar to that of entanglement ^IS] in the same 
model. 



2. Thermal state fidelity in strongly correlated systems 

The leading term of the fidelity, i.e. the thermal fidelity 
susceptibility is simply the specific heat if we choose the 
temperature as the driving parameter, and magnetic sus- 
ceptibility if we choose the magnetic field as the driving 
parameter. Despite of this there are still some works on 
fidelity in thermal phase transitions, especially for those 
Hamiltonians with non-commuting driving terms. 

The BCS superconductor: The effective Hamiltonian 
of the Ba rdeen-Cooper-Schrieffer (BCS) superconductor 
l21l| can be written as 

-ffecs = '^£kc\^Ck^ +'^Vkk'c\,^c^_f.,^c-kiCk,^. (257) 

k(T kk' 

where is the dispersion and Vkk' are coupling constant. 
The BCS Hamiltonian provides us with an example of 
a model with a mutually non-commuting Hamiltonian. 
The fidelity approach to its thermal phase transition from 
a normal state to a superconducting state was done by 




FIG. 28: (Color online) The gap A (UP) and the fidelity of 
SV = 10"^, (5T = (DOWN) as a function of the temperature 
T and the coupling constant V. The plot is given in rescaled 
quantities T kBT/{huJD), V DpV and A -> A/{hujD) 
(From Ref. [Hi). 

Paunkovic and Vieira [5^. Based on the mean-field ap- 
proach, the gap function and the thermal-state fidelity 
was obtained analytically. Fig. [28] shows both the gap 
and the fidelity as a function of the coupling and tempera- 
ture. Clearly, when the gap appears, the system becomes 
a superconducting state. A thermal phase transition oc- 
curs at the critical point. On the V — T plane, both V 
and T can be regarded as the driving parameter, there- 
fore, the fidelity between two thermal state separated by 
a small distance on the plane shows a sharp peak around 
the critical point. On the other hand, since the two terms 
at the right hand side of Eq. (|257p do not commute, the 
fidelity is not easy to be calculated. However, the au- 
thors found that distance obtained from Eq. (|187p and 
the fidelity have tha same qualitative bahavior (constant 
value 1 everywhere, except for the sudden drop along the 
line of phase transition) . 

Paunkovic and Vieira studied also the thermal phase 
transitions in the Stoner-Hubbard model whose Hamil- 
tonian reads 

Hs-a^'^Skcl^Cka + U'^n^nii. (258) 
where Sk = /2m is the dispersion, and U the on-site 
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interaction. Unlike the BCS superconductor, the phase 
transition they addressed in the Stoner-Hubbard model is 
due to the existence of Zeemann-like term for conserved 
quantities. In this case, the fidelity susceptibility is sim- 
ply the susceptibility of the conserved quantity 4§|. The 
fidelity shows a sharp peak around the transition point 
from an order phase at low temperature to a disorder 
phase at high temperature. 



VII. NUMERICAL METHODS FOR THE 
GROUND-STATE FIDELITY 



This section is for students who are not familiar with 
the numerical diagonalization and the density matrix 
renormalization group technique. 



A. Exact diagonalization 

Direct diagonalization: One of the main goals of quan- 
tum mechanics is to diagonalize Hamiltonian matraces. 
Except for a few cases, however, the Hamiltonian of most 
quantum many-body systems cannot be diagonalized ex- 
plicitly. Fortunately, advance in digital computers make 
it possible to diagonalize the Hamiltonian numerically. 
Up to now, many numerical methods rooting in com- 
puter science have been extensively developed, among 
which the most straightforward method is the exact di- 
agonalization. Though the system size that can be diag- 
onalized numerical is still not large, the numerical results 
actually are very instructive for studies of the quantum 
many-body systems. 

In the exact diagonalization, we need to express the 
Hamiltonian in a set of basis, and then diagonalize it 
numerically. As a warm-up example, let's consider a spin 
dimer system. The Hamiltonian of a spin dimer with the 
Heisenberg interaction reads 
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To see its matrix form, we use the eigenstates of op- 
erators as basis, i.e., {| tt), | Ti), | it), I ii)}- In this set 
of basis, the Hamiltonian can be written as 
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With some standard technique in Linear algebra, we can 
find that its ground state is a spin singlet state 
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[| n>-UT>], 
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with eigenvalue Eq 
states are 

*2 



-3, while three degenerate excited 



[| n> + UT>]. 



1 

ITT), *3 = Ui) 
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with higher eigenvalue £'1(2,3) — 1- 

Numerically, we can construct the basis by a set of 
integer which is usually a binary digit. For the above 
example, we can use 0, 1, 2, 3 (that are 00, 01, 10, 11 in 
binary system) to denote the basis | TT), I Ti), I TT), I ii), 
respectively. Then, we can define an array to save the 
Hamiltonian. There are some standard methods |212j | to 
diagonalize a matrix, such as Householder method and 
QR algorithm, or standar d lib raries, such as Linear Al- 
gebra Package(LAPACK) [Hll and Intel's Math Kernel 
Library (MKL) [21^ . 

Lanczos method: For a small size sample, such as a 
10 1/2-spin chain, the dimension of the Hamiltonian ma- 
trix is not very large. We are able to diagonalize the 
whole matrix. All eigenstates and eigenvalues can be ob- 
tained. While if the dimension of the matrix size is very 
large, say 100000, it is almost impossible to implement 
the traditional diagonalization methods. However, if we 
are only interested in the ground state of the Hamilto- 
nian, it is still possible for us to find the wavefunction 
via the famous Lanczos method. 

Here we would like to introduce the Lanczos method 
to diagonalize a large sparse matrix and how to calculate 
the fidelity and its susceptibility via the Lanczos method. 
We are not going to address the principle of the Lanczos 
method, instead we show that the perturbation nature of 
the method facilitates calculation of the fidelity suscep- 
tibility numerically. Now we briefly introduce the basic 
steps of the Lanczos method. 

The first step in the Lanczos method is to transform a 
large spare matrix H to a tri-diagonal matrix T as 



/ ai Pi 
Pi a2 



T„ 
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where m is the cut-off number. The cut-off number is de- 
termined by the precession. For a real symmetric matrix 
_ff , we need an orthogonal matrix to do such a transfor- 
mation, i.e. 



V 



which satisfy 



and 
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Therefore, starting from an initial (random) vector v, we 
have the foUowing relations 



Vl 

/3lV2 

/33V4 



V, 

Hvi - (vi|i/|vi)vi, 

HV2 - (V2|i?|v2)v2 - (Vi|i?|v2)vi, 
UWi - (V3|i?|v3)v3 - (V2|i?|v3)v2, 



/3m V, 



m+1 



-(Vm-l|i/|v,„)v„i_i. 



(268) 



From these relations, we can find that, using the orthog- 
onal condition (|266p . 



ai = (vi|iJ|vi), 




|r,) = (i/-a,)|v,) 


3i-i|vi-i 






v,+i) = |r,)/ft. 




= (v,j+i|i?|v,j+i). 




i — 1, 2, m — 1. 





(269) 

The above steps are the famous Lanczos iteration, which 
was proposed by Cornelius Lanczos in 1950. 

After the matrix T is obtained, one can easily calculate 
its eigenvalues Ei and their corresponding eigenvectors 
|ui). This procedure is simple because T is already tri- 
diagonal. It can be proved that the lowest eigenvalue of T 
is the ground-state energy of H. Then the ground-state 
wavefunction (Ritz eigenvector) can be calculated as 



\^a)^V\uo), 



(270) 



where jug) is the ground state of T. 

To calculate the fidelity between two ground states in 
parameter space, the Lanczos method is very promising 
because it is actually a numerically perturbative method. 
To be concrete, if i7(A)|^'o(A)) = £;o|^'o(A)), then for the 
Hamiltonian 

H{X + SX) = H{X) +SXH1, (271) 

we can use |5'o('^)) as an initial state, then 

ai =So + '5A(*o(A)|ff/|*o(A)), 

which is simply the ground-state energy up to the first- 
order perturbation. To the second order, one can find 
that 

ai = £;o + a(*o(A)|-ff/|*o(A)), 

ft = (5A^(vI/o(A)|F||*o(A)) - (vI/o(A)|i/,|vI/o(A))2, 

-^(*o(A)| Ai//i/(A + 6X)AHi\^oW), 
Pi 

(272) 



a2 



TABLE IL The fidelity susceptibility of the one-dimensional 
transverse-field Ising model obtained by the exact diagonal- 
ization (middle) and exact analytical results (right). Small 
difference might be caused by numerical derivation. Here 
N = 20. 



h 


FS(ED with Sh = 0.005) 


FS of Eq. (fml) 


0.2 


1.30206 


1.30208 


0.3 


1.37360 


1.37362 


0.4 


1.48807 


1.48809 


0.5 


1.66672 


1.66675 


0.6 


1.95535 


1.95539 


0.7 


2.48564 


2.48568 


0.8 


3.81093 


3.81096 


0.9 


7.96763 


7.96833 


1.0 


11.87198 


11.87500 


1.1 


5.83999 


5.84015 


1.2 


2.28066 


2.28060 


1.3 


1.13390 


1.13388 


1.4 


0.67720 


0.67719 


1.5 


0.44735 


0.44735 


1.6 


0.31372 


0.31372 


1.7 


0.22904 


0.22904 



where AHi = Hj - (^'o(A)|i7/|^'o(A)). Then the T ma- 
trix, to the second order, becomes 



T2 = 




The eigenvalue of T2 is 

1 

2 



ai + a2 - \ {ai - + 4/3f 
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Therefore, the Lanczos method is not only a perturbative 
method, but also a variational method. The advantage 
is that we can use 1^*0 (A)) as the initial state to calculate 
|4'o(A -I- SX)). Once the latter is obtained, the fidelity 
(^o(A)lvI'o(A -I- SX)) can be easily calculated. 

Here, we take the one-dimensional transverse-field 
Ising model as an example. In section IIIIl we can find 
the fidelity susceptibility at the point takes the form 



XF 



E 

k>0 



where 



dX 



d9k 
dX 



sin k 



1 



2 1 -f /i2 - 2ft, cos fc ' 



(275) 



(276) 



From Eq. I275i we can calculate the fidelity susceptibility 
up to a very large systems. For a comparison, we simply 
show both the analytical results and numerical results for 
a 20-site system in Table. |lTl 
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FIG. 29: (Color online)The infinite-system algorithm of tlie 
DMRG teclmique: The system size grows under the renor- 
malization group transformation. 



gorithm is more efHcient for a system with open bound- 
ary conditions. So we limit our discussion to this case 
here. Readers interested in other details of the DMRG 
algorithm are recomm ende d to refer to the review pa- 
per by U. SchoUwock (217j | and the book by Peschel et 
al |216l | . There are also some new development in the 
time-dependent DMRG method[2T8l|, and DMRG algo- 
rith m fro m the perspective of quantum information the- 
ory {219^ . 

The infinite-system algorithm (Fig. I29p : 

Step 1: Start with a small system block of I sites. Con- 
struct a set of real-space basis (for instance, {| t), | i)} for 
a single spin or {| tt), | Ti), I it), I ii)} for two spins). 
Then under the basis, construct the Hamiltonian ma- 
trix and matrices of operators responsible for inter- 
actions, say , at the rightmost site. 

Step 2: Construct the environment block of the same 
size (range from Z -I- 1 to 2/), including the Hamiltonian 
matrix and interaction-operator matrices, say Of, 
at the leftmost site, say O/, in a similar way. 

Step 3: Build the superblock by connecting the system 
and environment blocks. The superblock Hamiltonian 
can be constructed from 



+ O? ® Of 



B. Density matrix renormalization group 

The density matrix renormalization group (DMRG) 
technique is a very successful numerical method for one- 
dimensional strongly correlate d sys tems. The algorithm 



(277) 



was invented by S. R. White |215l | in 1992. Tzeng and 
Yang [ssj firstly applied the DMRG technique to calcu- 
late the ground-state fidelity. In this section, we briefly 
explain the basic ideas of the DMRG algorithm and its 
application to the fidelity. 

The main idea of the DMRG algorithm is a smart re- 
duction of the number of effective degrees of freedom and 
a variational and perturbative search of the ground state 
within the reduced space. This idea is very important in 
computational physics. For a quantum many-body sys- 
tem, the dimension of the Hilbert space grows exponen- 
tially with the system size. For example, a L-site spin- 
1/2 system, the dimension of its Hilbert space is 2^. The 
growth in the Hilbert-space dimension quickly exhausts 
computing resources. This fact discourage us from doing 
exact diagonalization for larger systems. Fortunately, the 
DMRG algorithm is able to capture the most relevant de- 
grees of freedom, hence allows us to reduce the effective 
dimension of the Hilbert space significantly. 

The DMRG algorithm consists of two fundamental sec- 
tions: infinite-system algorithm and finite-system algo- 
rithm. In both algorithms, a system of L sites usually is 
divided into four blocks, i.e, L'^-site system block S, two 
intermediate sites, and the rest environmental block E of 
— L — — 2 sites. The infinite-system algorithm 
aims to grow the system to the size we want to study, 
and the finite-system algorithm to reduce the numerical 
error based on the variational principle. The DMRG al- 



where — Of. ® Of is the interaction Hamiltonian be- 
tween the Ith. site of the system block and {I -\- l)th of the 
environment block. If there are more interaction terms, 
they should be included in too. Diagonalize the su- 
perblock Hamiltonian H to obtain the ground state |^'o) 
by the Lanczos method or Davidson method [220j |. Cal- 
culate the reduced-density matrices of the system block 
and the environment block by 



trB|*o)(1'o|,P^ = trs|*o)(*o| 



(278) 



then diagonalize by the traditional diagonalization 
method for dense matrix, such as Houscholdcr-QR 
method, 



w,- V,- 



(279) 



where Wj is in a decreasing order. Form a new set of 
basis for the system block by M eigenstates of with 
the largest eigenvalues, and construct the transformation 

matrix 



VL = {Vl, V2, • • • vm}. 



(280) 



Here M is chosen based on the desired precision. Trans- 
form the Hamiltonian of the system block and the inter- 
action operators at the rightmost site into the new basis 



of = n^ofn. 



(281) 
(282) 



The environment block Hamiltonian and Oi are trans- 
formed in a similar way. 
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Step 4: Connect a new site to the rightmost site of the 
system block, and a new site to the leftmost site of the 
environment block. The new Hamiltonian is 

= ®I + 1^ + H' , (283) 

= I (E) H'^ + H'^ + H' , (284) 

where H'^ is the single-site Hamiltonian. If the total size 
of the system block, environment block, and two sites is 
small than the target size, then goto the step 3, otherwise, 
jump to the following finite-system algorithm. 

The finite-system algorithm (Fig. I30|) : 

Once the size of the superblock in the infinite-system 
algorithm reaches the target size we want to study, one 
can calculate the ground-state properties. However, the 
results usually are not satisfactory because the error is 
still very large. We need to use the finite-system algo- 
rithm to reduce the error. The finite-system algorithm 
is similar to the infinite-system algorithm. Instead of 
growing both blocks in the infinite-system algorithm, the 
growth of one block is accompanied by shrinkage of the 
other block in the finite-system algorithm. The informa- 
tion of the block shrunk should be retrieved. For this 
purpose, one need to store the block information, includ- 
ing matrices of the Hamiltonian and interaction opera- 
tors in the reduced space, obtained in the infinite-system 
algorithm. 

Step 1: Followed from the infinite-system algorithm, 
now we have a system block of L/2 — 1 sites, two in- 
dependent sites, and an environment block of L/2 — 1 
sites. Construct the superblock Hamiltonian, and com- 
pute the ground state by the Lanczos method or David- 
son method. Differ from the infinite-system algorithm, 
here we do the reduced basis transformations for the sys- 
tem block only. Then the new system block has L/2 
sites. 

Step 2: Continue to grow the system block in a similar 
way, while using the stored information of the environ- 
ment block, until the system block reaches the maximum 
size. 

Step 3: Once the environment block is minimized, then 
grow the environment block at the expense of the system 
block in the same way, until the environment is maxi- 
mized. 

Step 4: Grow the system block, and shrink the envi- 
ronment block to L/2 — 1. 

A complete shrinkage and growth sequence for both 
blocks (From step 1-4) is called a sweep (Fig. [30l). Usu- 
ally more sweeps, higher precision the final results have. 
Once the desired precision is reached, the ground state 
can be expressed in reduced space. 

In the DMRG algorithm, the reduced basis obtained 
by the numerical renormalization group strongly depend 
on the parameter of the Hamiltonian. So we cannot com- 
pare two ground states at distinct points in the parameter 
space directly. In order to calculate the fidelity between 
the two ground states, one has to find a transformation 
between the two sets of reduced basis. Such a transfor- 



mation, as proposed in Ref. [70|, can be established in 
the final sweep of the finite-system algorithm. 

The ground-state wavefunctions of the Hamiltonian 
Hq + XHj at two points Ai and A2, in their own reduced 
space, can be expressed as 

|vI'o(Ai)> - 'i>.™n,(A)|<^f)|m)|n)|<^f), (285) 

i.'rn,n.j 

|*o(A2)) - J2 ^^mnJ{XM)\m)\n)\^f), (286) 

where \'pf){\'pf)) is the reduced basis of the sys- 
tem(environment) block, and |m), |n) are the basis of 
the middle two sites. Since |m), |n) are the basis of local 
sites, we have 

(m|m) = 5mm, {n\n) = 6nn- (287) 
Then the fidelity between two ground states becomes 

(^o(Ai)|*o(A2)) (288) 

where T^f, = {(pf\'(pf,) (Tfj, = {ipf\lpf,)) defines the trans- 
formation matrix between the two sets of reduced basis 
of the system (environment) block at Ai and A2. 

Now we focus on how to obtain the of a L'^-site 
system block, the of the environment block can be 
obtained in a similar way. 

Step 1: During the final sweep of the finite-system al- 
gorithm of the two Hamiltonians H{Xi) and H{X2), if 
both system blocks are minimized (to a single site) , their 
basis are defined in the real space and not reduced, then 
the transformation matrix between the two system blocks 
is simply the unity matrix, i.e., 

T^(l = 1) = /. (289) 

Step 2: Suppose the transformation matrix of the two 
system blocks up to I sites (including Z = 1) is obtained, 
i.e. 

T^{l) = \Tpf,{l)){^f{l)\, (290) 

the basis of the system blocks together and the new site 
before the renormalization group transformation are 

\ipf{l))®\mil + l)), (291) 
\^f,il))^\mil + l)), (292) 

respectively. The transformation matrix between the two 
sets of basis becomes 

T'^{l)(g> 1(1 + 1). (293) 

After the RG transformation, one can find the transfor- 
mation matrix becomes 

T'^{l + l) = n^[T'^{l)(g)I{l + l)]n, (294) 
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FIG. 30: (Color online)The finite-system algorithm of the 
DMRG technique: The system block grows under the renor- 
malization group transformation, under the environment 
block shrinks with the retrieved representation; and vice 
versa. 



ment in the ground state of quantum many-body sys- 
tems, the fidelity approach to quantum phase transitions 
have shed new light on the critical phenomena. Moreover, 
unlike the entanglement, which is somehow still mysteri- 
ous in quantum many-body system, the fidelity has, from 
our point of view, a clearer physical picture. Especially, 
its leading term manifests scaling and critical behaviors 
around the phase transition point. Therefore, the fidelity 
is really a new optional method to investigate quantum 
phase transitions, especially for those cases that we know 
nothing about order parameters. 

However, there are still some remaining problems. 
Firstly, the validity of the fidelity in the Beresinskii- 
Kosterlitz-Thouless phase transitions is still controver- 
sial. Secondly, the deep reason that the fidelity can sig- 
nal the topological phase transitions remains unknown. 
Thirdly, one of most difficulties of the fidelity in quan- 
tum many-body systems is that it is really difficult to 
find the exact ground state, except for a few cases. Fi- 
nally, it is a challenging problem to measure the fidelity 
in experiment on scalable systems. 

Finally, since the field is still quickly developing, we 
hope again that this introductory review can offer some 
rough essays first, then to arouse other people's better or 
more mature ideas. 



which gives a recursion relation of the transformation 
matrix in the RG transformation. 

Step 3: repeat the step 2 until both the two system 
block grow upto sites. 

Clearly, the transformation matrix can be obtained 
if we start from the minimized environment block. Fi- 
nally, the fidelity of Eq. (|288[) can be calculated. The 
fidelity susceptibility can be also computed by the nu- 
merical differentiation. 



VIII. SUMMARY AND OUTLOOK 

As mentioned in the introductory section, quantum in- 
formation theory provides us opportunities to investigate 
quantum phenomena from new angles. It is fair to say, 
besides huges of works on the role quantum entangle- 
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